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INVESTIGATION  OF  UNDERGROUND  EXPLOSIONS  WITH  MODEL  TESTS: 
MEASUREMENTS  IN  THE  TUBE 


SUMMARY 


Model  tests  in  three  scales  were  undertaken  to  examine 
if  simple  scaling  laws  can  be  used  to  investigate 
underground  explosions.  Charges  were  detonated  in 
steel  tubes,  and  pressures  and  arrival  times  we  re 
measured  in  the  tube  and  on  a platform.  This  report 
describes  the  results  from  the  measurements  in  the 
tube,  A systematic  deviation  from  the  scaling  lav/s 
is  observed.  The  dependence  of  this  deviation  on  tube 
diameter,  charge  weight  and  wall  roughness  is  examined. 

The  results  from  the  measurements  on  the  platform  are 
given  in  Ref.  (1) . A detailed  description  of  the  experi- 
mental 'fidget  is  given  in  Ref  (2)  and  of  the  data  pro- 
cessing in  Ref  (3) . 


INTRODUCTION 


This  report  describes  the  results  from  tests  performed  to 
examine  if  simple  scaling  laws  can  be  verified  for  the  pre- 
diction of  the  blast  wave  following  an  accidental  exolosion 
in  an  underground  ammunition  storage  site.  Such  scaling  laws 
have  been  verified  for  explosions  in  free  air,  and  it  has 
been  proposed  that  simple  scaling  laws  also  may  be  used  in 
the  underground  case  (4),  (5),  (6). 

There  are,  however,  many  effects  which  obviously  do  not  scale 
according  to  these  laws.  These  effects  may  be  of  such  a low 
order  of  magnitude  that  for  practical  purposes  they  are  negli- 
gible. This  has,  however,  not  yet  been  sufficiently  verified. 

The  objective  of  our  tests  was  to  compare  the  shock  parameters, 
front  pressure  and  front  velocity,  when  TNT-charges  of  diffe- 
rent weights  were  detonated  in  models  in  different  scales , 
and  to  see  if  there  was  any  deviation  of  importance  from  the 
scaling  laws.  Models  of  a simplified  underground  storage 
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site  were  therefore  built  in  three  scales.  The  models 
consist  of  a tube,  simulating  a tunnel,  and  a platform, 
simulating  a flat  terrain.  The  scales  are  1:100,  2:100 
and  3:100  of  a natural  size.  TNT-charges  were  detonated 
and  the  pressure  time  history  of  the  resulting  blast 
was  recorded  at  measurement  stations  in  the  tube  and  on 
the  platform. 

The  results  from  the  measurements  inside  the  tube  are 
described  in  the  present  report,  while  the  results  from 
the  measurements  on  the  platform  are  described  in  Ref  (1). 
The  experimental  equipment  is  described  in  Ref  (2)  together 
with  the  methods  for  recording  the  data  on  analog  tapes. 

The  data  are  digitalized  and  are  processed  on  a CYBER  74 
computer.  This  is  described  in  Ref  (3) . 

Some  of  the  results  from  the  tests  have  been  published  in 
a preliminary  report  (4) . 


2 SCALING  LAWS 


In  Figure  2.1  is  shown  a sketch  of  the  tube  with  the  charge. 
The  tube  has  diameter  D.  The  cylindrical  charge  has  a dia- 
meter vD  and  length  d.  The  charge  weight  Q is  then  propor- 
tional to  dD^.  The  pressure  in  the  undisturbed  air  ahead 
of  the  shock  front  is  Pq,  and  the  pressure  immediately  behind 
the  front  is  P2.  If  only  these  parameters  are  used  to  de- 
scribe the  problem  it  follows  from  dimensionless  analysis  that 
the  front  overpressure 


in  a distance  L from  the  charge  center  can  be  expressed  as 
a function  of  the  form 
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(2.2) 


If  we  introduce  the  charge  weight  Q in  this  equation,  we 
obtain 


'ld2  q 
Q ' d3 


(2.3) 


It  will  be  convenient  to  introduce  a parameter  q,  called 
charge  group,  defined  by  the  equation 


q 


(2.4) 


It  is  also  convenient  to  introduce  a dimensionless  parameter 
for  distance,  namely 


n 


(2.5) 


n measures  the  distance  of  the  shock  front  from  the  charge 
center  in  number  of  tube  diameters. 


With  n and  q introduced  the  scaling  law  takes  the  form 


P = 


q 


(2.6) 


It  states  that  for  each  tube  diameter  D the  front  pressure 
will  the  same  at  trie  same  position,  when  measured  in  numbe 
of  tube  diameters , and  when  the  charges  detonated  belong  to 
the  same  charge  group,  i e have  the  same  value  of  q.  We  shal 
call  this  scaling  with  respect  to  dimension. 

Since  q is  proportional  to  ~ it  describes  the  geometrical 
conditions  in  the  vicinity  of  the  charge * In  a position  not 
too  close  to  the  charge,  it  is  a reasonable  assumption  that 
this  term  in  the  argument  of  f may  be  neglected , leaving  the 
scaling  law  on  the  simple  form 
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P = f (Lg) 

where 


(2.7) 


(2.8) 


is  called  scaled  distance. 

For  the  same  tube  diameter  the  front  pressure  will  be  the 
same  for  different  charge  groups  when  it  is  measured  in 
the  same  scaled  distance  Ls , i e when  ^ is  held  fixed. 

This  is  called  self  similarity. 

By  similar  arguments  as  above  the  arrival  time  t of  the 
shock  front  in  a position  L is  given  by  a functional  relation 
of  the  form 

ts  = 9(Ls'  5)  (2.9) 

where  we  have  introduced  the  scaled  time 


At  a sufficiently  large  distance  from  the  charge,  this 
scaling  law  is  assumed  to  reduce  to  the  simple  relation 

ts  - g(Ls)  (2.11) 

The  scaling  laws  (2.7)  and  (2.11)  were  proposed  by  Erikson 
(5)  in  1964  for  the  propagation  of  blast  in  a tube. 

The  validity  of  these  scaling  lav/s  will  be  examined  in 
this  report. 
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3 DESCRIPTION  OF  THE  EXPERIMENTS 

3 . 1 Description  of  the  models 

The  models  consist  of  two  parts.  One  is  the  tube  which 
simulates  the  underground  tunnel.  The  other  is  the  platform 
section  consisting  of  a vertical  steel  plate  above  the  tube 
outlet  and  a horizontal  platform,  simulating  a steep  mountain 
side  and  a flat  terrain  respectively. 

The  tubes  are  designed  to  provide  two  different  model  con- 
figurations . 

The  first  configuration  (CFG1)  consists  in  principle  of  a 
straight  tube  with  a detonation  section  in  the  middle  as 
shown  in  figure  3.1. 

The  second  configuration  (CFG2)  in  principle  consists  of  a 
straight  tube  connected  to  a detonation  chamber  in  the  one 
end.  The  other  end  leads  to  the  platform  as  shown  in 
frgure  3.2.  The  cross  section  area  of  the  detonation  chamber 
is  twice  the  cross  section  area  of  the  tube. 

The  tunnel  diameters  in  both  configurations  are  0.05  m, 

0.10  m and  0.15  ra  respectively  for  the  scales  1:100,  2:100 
and  3:100. 

Both  models  were  made  to  allow  for  the  mounting  of  pressure 
transducers  along  the  tubes,  one  for  each  four  diameter  length. 
Pressure  transducers  were  also  mounted  on  the  platform  along 
the  projection  of  the  extended  tube  axis  with  one  per  each 
0.5  m approximately.  See  figures  3.3  and  3.4. 

Both  configurations  allow  for  a large-scale  roughness  in 
the  tubes  to  take  three  different  values,  rL  = 0 , rL  = 0.03 
and  mp  = 0.06.  The  large-scale  roughness  rL  is  defined  by 
the  relation 
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(3.1) 


D is  the  tube  diameter  and  2 AD  is  the  height  of  the  roughness 
element  as  shown  in  figure  3.5. 

The  charges  are  pressed  TNT  charges.  They  are  of  cylindrical 
shape  and  for  both  configurations  the  diameters  are 
the  same  fraction  v of  the  tube  diameter  D in  the 
three  scales.  v = 0.38  for  CFG1  and  v = 0.29  for  CFG 2 . 

There  is,  however,  some  difference  in  the  mounting  of  the 
igniters,  as  is  described  in  further  detail  in  Ref  (2). 

The  measurements  and  recordings  of  the  data  are  automatic. 
This  is  further  described  in  Ref  (2)  and  (3) . 


3 . 2 Rankine-Hugoniot  equation 

We  have  measured  the  arrival  time  of  the  shock  front  and 
the  front  pressure  at  different  measurement  stations.  From 
the  arrival  time  measurements  the  shock  front  velocity  U 
can  be  calculated.  The  shock  velocity  U and  the  front  over- 
pressure P are  related  to  each  other  by  the  Rankine-Hugoniot 
equation 


P 


(1  + y2) 


1 


(3.2) 


aQ  is  the  speed  of  sound  in  the  undisturbed  air  ahead  of 
the  shock  front  and 


where  y is  the  adiabatic  constant.  This  relation  applies 
to  ar.  ideal  gas.  The  value,  y = 1.4  for  air  is  used  in  our 
calculations . 
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By  this  formula  then,  the  pressure  and  arrival  time  measure- 
ments represent  a double  set  of  data. 


3 . 3 Experimental  series 

The  main  series  of  experiments  consisted  of  approximately 
three  shots  in  each  of  six  charge  groups  in  the  three  scales 
and  for  the  three  values  of  wall  roughness  and  the  two  con- 
figurations. The  charge  groups  and  the  corresponding  charge 
weights  in  the  three  scales  are  listed  in  Table  3.1  together 
with  the  corresponding  loading  densities. 

The  shots  with  the  largest  charges  in  scale  3:100,  CFG2, 
were  for  different  reasons  only  performed  in  the  case  of 
rL  “ 0 . 

It  was  of  interest  to  perform  some  shots  with  charge  groups 
corresponding  to  smaller  loading  densities.  This  was  done 
in  sea] e 3:100  with  r ^ = 0 for  both  configurations.  The 
charge  weights  used  are  listed  in  Table  3.2.  Three  shots 
were  fired  for  each  of  the  charge  weights  in  each  configu- 
ration . 


3 • 4 Short  description  of  the  propagation  of  the  bl a st 

The  cnarges  are  placed  coaxially  with  the  tubes.  When  the 
charge  is  detonated  a blast  wave  will  propagate  through  the 
tube  and  onto  the  platform.  The  pressure-time  history  and 
the  arrival  time  of  the  shock  front  are  measured  and  recorded 
at  every  measurement  station  in  the  tube  and  on  the  platform. 

Figure  3.6  shows  as  an  example  the  front  pressure  in  the  13 
measurement  stations  in  the  tube  for  three  equal  shots  in 
scale  1:100,  CFG1.  We  notice  the  low  pressure  close  to  the 
cnarge.  The  shock  front  needs  some  time  to  build  up.  This 
is  partially  due  to  the  fact  that  the  energy  is  not  instant- 
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aneously  released  and  that  the  initial  expansion  of  the  shock 
is  approximately  spherica] . First  after  a set  of  reflections 
from  the  walls  a front  is  built  up.  The  distance  at  which 
the  front  is  established  varies  with  charge  weight.  At  the 
fourth  measurement  station,  16  diameters  from  the  charge 
center,  it  is  always  fully  established.  We  shall  denote 
these  first  16  diameters  the  build-up  region  and  treat  it 
separately . 

From  position  n = 16 , where  the  shock  is  established,  the 
shock  wave  may  be  described  as  one-dimensional  (5)  (6). 

In  this  region  we  shall  study  the  attenuation  of  the  front 
pressure  as  it  propagates  through  the  tube. 

We  see  from  Figure  3.6  that  the  pressure  does  not  decay 
smoothly  with  distance.  This  does  not  only  result  from 
uncertainties  in  the  pressure  measurements.  It  is  also 
due  to  physical  effects  such  as  reflections  from  the  tube 
walls  and  to  rarefaction  v; aves  overtaking  the  shock  front. 

As  a xesu.lt  of  this  complicated  system  of  reflections  and 
rarefaction  waves  we  get  the  relatively  large  scatter  of 
the  measurements  from  equal  shots,  as  seen  in  the  figure. 

This  scatter  is  observed  in  all  the  measurements.  We  shall 
call  the  region  between  position  n = 16  and  the  tube  outlet 
for  the  attenuation  region,  since  the  attenuation  of  the 
shock  is  the  dominating  feature. 

We  shall  not  include  the  measurement  station  in  the  tube 
outlet  in  this  region,  that  is  n = 52  for  CFG1  and  n = 60 
for  CFG2.  When  the  shock  waves  passes  through  the  outlet, 
it  changes  from  one-dimensional  to  three-dimensional.  This 
region  will  be  studied  separately  in  Ref  (1) , where  also 
the  attenuation  of  the  shock  on  the  platform  will  be  examined. 
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3 . 5 Systematic  errors 

It  is  not  easy  to  estimate  the  magnitude  of  the  uncertainties 
of  the  measurements,  but  it  is  reasonable  to  believe  that 
it  is  in  the  order  of  10%  for  the  pressure  measurements  and 
5%  for  the  arrival  time  measurements.  This  is  discussed 
in  Ref  (3)  . 


Two  systematic  errors  related  to  the  measurements  should  be 
mentioned. 


One  is  that  two  different  methods  have  been  used  in  recording 
the  front  pressures,  and  the  other  is  the  effect  of  the 
igniters.  They  will  be  shortly  described, 

3.5.1  Differences  in  recording  of  front  pressures  and  arrival 
times  between  the  scales 

The  pressure  time  history  of  the  blast  at  each  measurement 
station  was  recorded  on  analog  tapes  (2) . For  the  experi- 
ments in  scale  1:100  these  were  played  out,  and  the  values 
of  the  front  pressure  were  found  by  direct  measurements  on 
the  resulting  diagrams.  For  the  experiments  in  scales  2:100 
and  3:100  the  results  were  transmitted  to  a digital  tape, 
and  the  front  pressures  were  found  by  an  automatic  routine 
(3)  . 

For  a number  of  shots  both  methods  were  used,  and  it  was 
found  that  the  automatic  method  usually  gave  somewhat  lower 
values  for  the  front  pressures.  The  difference  was  of  magni- 
tude 2 bar.  As  a consequence,  the  recorded  front  pressures 
in  scale  1:100  will  have  a tendency  to  be  a little  too  high, 
compared  with  the  other  two  scales, 

A similar  comparison  showed  that  the  automatic  method  gave 
arrival  times  of  the  shock  about  0.01  ms  larger  than  the 
manual  method.  This  results  in  too  small  arrival  times  in 
scale  1:100  compared  with  the  two  larger  scales.  This 
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difference  will  be  of  importance  at  the  measurement  stations 
relatively  close  to  the  charge. 

3.5.2  TNT-equivalence  of  the  electric  igniters 

We  have  used  microsecond  electric  igniters  P no  8 from 
Nitro  Nobel  AB,  Sweden.  They  are  not  of  standard  type. 

For  the  smallest  charges  the  TNT-equivalence  of  the  igniters 
may  be  of  importance.  Some  shots  were  therefore  made  to 
determine  this. 

The  shots  were  fired  in  scale  1:100,  with  rL  = 0 and  in 
CFG2 . Front  pressures  and  arrival  times  were  recorded  at 
the  13  measurement  stations  in  the  tube.  From  the  arrival 
times  the  front  pressures  were  calculated  by  the  Rankine- 
Hugor.oit  equation.  The  front  pressures  determined  from  the 
arrival  times  were  in  good  accordance  with  the  measured 
front  pressures. 

The  front  pressures  resulting  from  a charge  consisting  of 
10  g TNT  and  one  electric  igniter  were  compared  with  the 
front  pressures  from  a charge  consisting  of  10  electric 
igniters.  In  the  same  manner  a charge  consisting  of  20  g 
TNT  and  one  ingiter  was  compared  with  a charge  consisting 
of  20  igniters.  There  were  fired  2 shots  with  10  g TNT, 

3 shots  with  10  igniters,  2 shots  with  20  g TNT  and  2 shots 
with  20  igniters.  The  results  are  listed  in  Table  3.3. 

If  P‘  denotes  the  pressures  resulting  from  TN'T-charges  and 
F"  the  pressures  from  the  igniters,  it  was  found  that  ~ 
was  approximately  constant  at  the  different  measurement 
stations  in  the  tube.  For  the  shots  with  10  g TNT  and 
10  igniters  it  was  found  to  be  —7-  = 0.7  8 for  the  measured 

p II  P 

pressures  and  —7  --  0.79  for  the  calculated  pressures.  For 
the  shots  with  20  g TNT  and  20  igniters  it  was  £-7-  = 
for  both  the  measured  and  the  calculated  pressures. 


0.81 
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For  the  order  of  magnitude  analysis  in  which  we  are  inter- 
ested we  shall  assume  that  the  pressures  measured  are  pro- 
portional to  the  charge  weight.  If  x denotes  the  TNT- 
equivalence  of  the  igniters  we  obtain 


P"  = 10  x 

P'  10  g + x 


(3.4) 


for  the  case  of  10  g TNT  and  10  igniters,  and 


P^  = 20  x 

P'  20  g + x 

for  the  case  of  20  g TNT  and  20  igniters, 
x = 1.15  g 


(3.5) 


This  gives 

(3.6) 


in  the  first  case  and 

x = 1 . 0 7 g (3.7) 

in  the  second  case. 


This  shows  that  the  effect  of  the  igniters  is  of  some 
importance  for  the  10  g charge  and  the  20  g charge.  The 
10  g charge  is  effectively  an  11  g charge  and  the  20  g charge 
is  effectively  a 21  g charge,  representing  an  increase  of 
10  and  5 percent,  respectively.  Therefore,  when  the  results 
from  these  shots  are  compared  to  the  shots  with  the  same 
charge  groups  in  the  two  larger  scales , the  pressures  in 
scale  1:100  are  somewhat  too  high. 

This  effect  will  be  of  no  importance  for  larger  charge  groups 
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3.6  Short  review  of  the  results 


As  mentioned,  the  main  series  of  shots  consisted  of  3 shots 
for  each  combination  of  the  different  charge  groups,  scales, 
configurations  and  degrees  of  roughness.  This  gives  a total 
of  about  350  shots.  For  every  shot  the  arrival  time  of  the 
shock  front  of  the  blast  and  its  pressure  time  history  were 
recorded  at  each  measurement  station.  The  front  pressure 
and  maximum  pressure  were  specially  extracted  from  the 
pressure-time  history.  The  work  reported  here  is  mainly 
based  upon  the  study  of  front  pressure,  maximum  pressure 
and  time  of  arrival  of  the  shock  front  at  each  of  the  3 
measurement  stations  in  the  build-up  region  and  the  9 in 
the  attenuation  region.  The  conclusions  presented  in  the 
following  are  thus  based  upon  a total  amount  of  data  on 
about  350  * (3+9)  * 3 = 12  600. 

3.6.1  Deviation  with  respect  to  dimension 

Let  us  as  an  example  look  at  the  results  for  one  charge 

3 

group,  q = 160  kg/m  . Figure  3.7  shows  the  mean  values 
of  the  front  pressure  in  the  build-up  and  attenuation  region 
for  the  three  scales.  According  to  the  scaling  lav;,  eq  (2.7) 
this  should  be  three  overlapping  curves.  At  first  sight 
these  curves  may  seem  to  be  equal  when  compared  to  the  scatter 
of  the  pressures  from  equal  shots,  as  shown  in  Figure  3.6. 

As  a mean  value,  however,  the  pressure  increases  with  in- 
creasing scale. 

Figures  3.8,  a-d,  show  the  front  pressure  as.  s function  of  tube 
diameter  for  some  charge  groups  and  positions.  According 
to  the  scaling  law,  these  should  be  straight,  horizontal 
lines.  Although  we  can  observe  all  kinds  of  curve  forms, 
on  the  whole  the  pressure  increases  with  increasing  scale. 

The  question  is,  however,  whether  the  increase  is  so  strong 
compared  to  the  data  scatter  that  it  must  be  concluded  that 
the  scaling  law  is  not  valid  according  to  our  data.  To 
answer  this,  a systematic  statistical  examination  is  necessary, 
which  will  be  presented  in  the  next  chapter. 
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3.6.2  Deviation  from  self-similarity 

According  to  self -similarity  one  should  for  each  scale 
measure  the  same  front  pressure  when  the  scaled  distance 
Lg  = “ is  fixed.  Figures  3.9,  a-c , show  the  front  pressure  as 
a function  of  charge  group  for  some  values  of  Ls  in  different 
scales.  These  should  be  straight,  horizontal  lines.  Again, 
we  observe  all  kinds  of  curve  forms,  but  on  the  whole,  the 
pressure  decreases  with  increasing  charge  group.  It  thus 
seems  that  an  increase  in  distance  from  the  charge  is  not 
sufficiently  compensated  for  by  increasing  the  charge  weight 
with  the  same  factor. 

Again,  however,  a systematic  statistical  analysis  is  neces- 
sary to  decide  if  the  deviation  is  a real  one,  or  if  it  might 
be  explained  by  the  scatter  of  the  measurements . 

4 LINEAR  REGRESSION 

4 . 1 Deviation  from  scaling . Students  T-test 

The  discussion  in  the  preceeding  chapter  showed  the  necessity 
or  a statistical  analysis  of  the  data.  It  will  be  presented 

t 

here . 

In  chapter  2 we  formulated  the  scaling  laws 

P = f(Ls)  (2.7) 


and 


ts  = g(Ls)  (2.11) 

We  shall  examine  the  possibility  of  a deviation  from  these 
scaling  laws,  and  we  let  the  deviation  be  a linear  expression 
in  D,  q and  n.  That  is, 


we  assume 
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P = f (Ls)  + ai  D + a2  q + a3  n 

= f(Ls)  + ai  D + [ a2  + a3  Ls]  q (4.1) 

and 


tg  = g(Ls)  + ai'  D + a2  ' q + 0,3  ' n 

= g (Ls)  + ax'  D + [ a2 ' + 03'  Ls] q (4.2) 

where  ap , a2,  a2  and  otp  ' , a2 1 , a 3 1 are  parameters  to  be 
determined  by  the  data.  ap  and  otp  ' measure  the  deviation 
from  Eqs  (2.7)  and  (2.11)  with  respect  to  scale  D. 

[ a2  - 013  Ls  ] and  [ a2 1 + 03'  Ls  ] measure  the  deviation  from 
self-similarity . 

The  _inear  form  does  of  course  not  necessarily  represent 
the  best  fit  to  the  deviation,  but  a possible  non-zero 
deviation  will  be  revealed.  For  the  scaling  laws  not  to  be 
denied  it  is  a necessary  condition  that  the  a-parameters  turn 
out  to  be  equal  to  zero.  More  precisely,  they  must  be  so  close 
00  zero  that  within  the  scatter  of  the  data  they  can  be  said 
to  be  equal  to  zero.  This  is  assessed  by  the  Students  T-test. 

We  shall  test  the  following  null-hypotheses: 

H0:  ci  = 0 H0:  a2  = 0 H0:  a 3 = 0 

HqS  ax'  = 0 H0 : a 2 1 = 0 H0:  a3 ' = 0 

The  test  level  is  e,  i e,  the  probability  that  the  hypothesis 
is  correct  although  we  by  the  test-criterion  has  denied  it, 
shall  be  less  than  e.  The  test  criterion  is  that  H0  is 
denied  if 


T > C 


(4.3) 
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where  T is  Students  T- statistic  for  the  actual  a-parameter. 


and  3 is  the 


1 - 


-fractil  for  the  Student  distribution 
with  N-2  degrees  of  freedom.  N is  the  number  of  observa- 
tions included  in  each  test,  which  in  this  connection  can 
be  set  equal  to  infinity. 


The  values  of  C for  some  values  of  test  level  e are  listed 
in  Table  4.1. 


In  the  next  sections  we  shall  give  the  results  of  these 
tests  for  the  build-up  region  and  the  attenuation  region. 


4 . 2 Build-up  region 

The  arguments  used  to  establish  that  P and  ts  are  functions 
of  n and  q only  through  the  scaled  distance  L«,  = - , is  valid 

q 

m a region  not  too  close  to  the  charge.  The  self-similarity 
is  in  other  words  not  supposed  to  be  valid  in  the  build-up 
region.  Since  there  are  only  three  measurement  stations  in 
this  region,  and  because  of  the  great  scatter  of  the  measure- 
ments , specially  in  position  n — 4 , the  validity  of  self- 
similarity is  not  examined. 

4.2.1  Deviation  with  scale 

It  is  important  to  examine  if  there  is  established  a deviation 
of  the  front  pressure  with  scale  already  in  the  build-up 
region.  In  Table  4. 2.  a are  shown  the  values  of  a ]_  for  the 
two  configurations  and  the  three  degrees  of  wall  roughness 
together  with  the  corresponding  values  of  T.  T > C in  all 
cases  for  a test  level  e = 2%,  and  in  four  of  the  cases 
T > C even  for  a level  e = 0.2%.  This  shows  that  the  hypo- 
thesis H0:  a 2 = 0,  must  be  denied. 

al  -^s  positive.  This  means  that  the  front  pressure  increases 
with  increasing  scale.  The  variation  within  our  three  scales 
is  about  10-20  bar,  and  since  the  magnitude  of  the  pressures 
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in  this  region  of  the  tube  is  about  100  bar,  this  deviation 
is  appreciable. 

The  values  of  ag1  giving  the  deviation  of  the  scaled  arrival 
times  with  tube  diameter  are  shown  in  Table  4.2b.  The  corre- 
sponding values  of  T are  also  shown.  In  CFG2  a]_'  is  negative, 
and  H0:  ai*  = 0 must  be  denied  since  T > C on  a test  level 
e = 0.2%  in  two  cases  and  e = 1%  in  one  case.  For  CFG1  it 
seems  that  ag'  =0  is  possible. 

As  described  in  section  3.5.1  the  difference  between  the 
recording  methods  resulted  in  smaller  values  for  the  arrival 
time  in  scale  1:100  compared  with  the  two  larger  scales. 

It  seems  likely  that  this  fact  explains  why  the  tendency 
for  the  arrival  times  is  not  so  clear  as  for  the  pressures. 

Thus,  it  occurs  to  be  a sound  conclusion  , when  the  results 
for  the  pressures  and  the  arrival  times  are  taken  into  account, 
that  there  is  established  a deviation  from  the  scaling  laws 
already  in  the  build-up  region.  The  front  pressure,  measured 
at  the  same  scaled  distance,  increases  with  increasing  scale. 


4 . 3 Attenuation  region 

4 

4.3.1  Deviation  with  scale 

We  shall  first  look  at  the  deviation  with  tube  diameter 
starting  with  the  pressure  measurements.  The  values  of  ag 
and  T are  shown  in  Table  4.3a.  T > C for  both  configurations 
and  all  values  of  rL  at  the  test  level  e = 0.2.  The  hypo- 
thesis H0:  ot]_  - 0 must  thus  be  denied.  Since  al  is  positive 
in  all  cases,  one  must  conclude  that  also  in  the  attenuation 
region  we  will  on  the  average  observe  higher  pressures  with 
increasing  scale. 

When  we  look  at  the  arrival  times,  we  find,  just  as  we  did 
for  the  build-up  region,  that  the  results  are  not  quite  clear. 
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Table  4.3b  gives  the  values  of  a^'  and  T.  For  CFG  1 with 
rL  “ ® ai  •’-s  positive,  while  in  the  other  cases  it  is  negative 
In  the  three  cases  in  CFG  2,  HQ:  ' = 0 must  be  denied  on  a 

level  e = 0.2%,  while  for  CFG  1 with  rL  = 0.06  it  must  be 

denied  with  e = 1%,  and  for  rL  = 0.03  it  must  be  denied  with 
a test  level  e = 20%. 

In  summary,  for  the  time  measurements,  H0r  a3  ' = o must  be 
denied,  and  the  arrival  time  is  found  in  most  cases  to  decrease 
with  increasing  scale,  this  being  consistent  with  the  obser- 
vation that  the  front  pressures  increase  with  increasing  scale. 

One  must  conclude  that  also  in  the  attenuation  region  the 
front  pressures  are  larger  in  the  larger  scales.  The  method 
presented  here  does  not  reveal  how  the  deviation  in  pressure 
with  scale  varies  in  the  different  positions  in  the  tube. 

This  will  be  studied  in  chapter  4 in  connection  with  the 
curvefitting  procedure. 

4.3.2  Deviation  from  self-similarity 

Table  4.4a  gives  the  values  of  and  together  with  the 
respective  values  of  T.  On  a level  e = 2%  the  hypothesis, 

H0;  a-2  = 0,  must  be  denied  in  all  the  cases.  The  situation 
is  not  so  clear  for  a3,  but  when  the  values  of  T are  compared 
with  the  values  of  C in  table  3.1  it  is  reasonable  also  to 
deny  the  hypothesis  H0:  a3  = 0.  In  every  case  [a3  + a3  Ls ] 
is  different  from  zero.  Since  both  a3  and  a3  are  negative, 

[a2  + a3  Ls ] is  negative,  and  the  deviation  AP  from  self- 
similarity 

AP  = q [a 2 + a3  Ls]  (4.4) 

j 

is  negative.  That  is,  the  pressure  decreases  when  the  charge 
increases,  although  the  scaled  distance  Ls  is  constant.  We 
also  notice  that  the  deviation  from  self similarity  is  larger 
for  larger  values  of  Ls. 
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The  time  measurements  give  similar  results.  The  values  of 
a2'  and  a3 ' together  with  the  corresponding  values  of  T are 
shown  in  table  4.4b.  Except  for  one  case,  H0:  a2  ' =0  and 
H0:  a 3 ' = 0 must  be  denied.  Although  a2 ' is  negative,  a3  ' 
is  positive  and  sufficiently  large  so  that  the  deviation  of 
scaled  time  from  self -similarity , 

Ats  = q (a2 ' + a3 ' Ls)  (4.5) 

is  positive  and  increasing  with  increasing  charge  weight, 
except  for  the  smallest  values  of  I,s.  This  positive  deviation 
in  scaled  arrival  time  is  equivalent  to  a negative  deviation 
in  front  pressure. 


4 . 4 Summary 

It  has  been  shown  in  this  chapter  that  the  scaling  laws, 
eqs(2.7)  and  (2.11)  neither  in  the  build-up  region  nor  in  the 
attenuation  region  are  confirmed  by  our  measurements.  The 
deviation  is  not  so  small  that  it  can  be  explained  by  the 
scatter  of  the  measurements. 

It  is  desirable  then,  to  describe  the  deviation  in  more  detail, 
trying  to  explain  the  reasons  for  the  deviation. 


5 CURVE  FITTING 


5.1  Introduction 


In  chapter  2 we  stated  the  scaling  law 

P = f (Ls)  (2.8) 

for  the  front  pressure,  preliminarily  neglecting  the  functio- 
nal form  f. 


21 


It  has  been  proposed  to  use  (5) 

P = f (Ls)  = Cg  LS"C2  (5.1) 

A curve-fitting  procedure  has  now  been  used  to  determine 
values  for  the  parameters  Cg  and  C2  from  the  experimental 
data  in  the  attenuation  region.  Since  we  have  observed  that 
the  scaling  law,  eq  (2.7)  is  not  satisfied  by  our  data,  we 
must  expect  Cg  and  C2  to  be  functions  of  both  scale  and  charge 
group.  We  also,  of  course,  expect  them  to  be  dependent  on 
wall  roughness  and  configuration. 

The  curve  fitting  is  performed  by  the  method  of  least  squares. 
This  is  done  separately  for  shots  belonging  to  the  same  con- 
figuration, wall  roughness,  scale  and  charge  group.  On  the 
average  then,  the  curve  fitting  is  based  on  the  measurements 
in  the  attenuation  region  for  three  equal  shots,  giving  27 
corresponding  values  of  pressure  and  position. 

The  results  are  listed  in  tables  5.1a  - 5. If.  They  show  the 
values  of  Cg  and  C2  for  each  configuration,  roughness,  charge 
group  and  scale,  together  with  the  correlation  coefficient,  R. 

We  first  notice  the  good  correlation  obtained,  and  that  it  is 
about  equally  good  for  each  configuration,  scale,  roughness 
and  charge  group.  For  example,  eq  (5.1)  does  not  seem  to  fit 
better  to  the  tubes  with  rL  --  0 than  to  those  with  rg,  = 0.03 
or  rL  = 0.06. 

We  notice  that  Cg  and  C2  are  functions  of  q and  D.,  and  dependent 
on  rg,  and  configuration's  they  should  be  according  to  the 
analysis  in  chapter  4. 

The  curve  fitting  procedure  will  now  be  used  to  study  the 
deviation  from  scaling  more  closely. 
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Since  we  have  divided  the  tube  into  two  regions,  the  build-up 
region  and  the  attenuation  region,  it  is  convenient  to  express 
eq  (5,1)  in  a slightly  different  form,  namely 


P 


o 


(5.2) 


PG  is  the  front  over-pressure  in  position  n = 16 , that  is  at 
the  beginning  of  the  attenuation  region.  The  pressure  at 
this  position  is  determined  by  what  has  happened  in  the  build- 
up region.  The  relation  between  PG  and  C] , C2  is 


o 


= Cn 


SL_ 

16 


C2 


(5.3) 


The  values  of  PQ  are  a] so  listed  in  the  tables  5.1. 


Since 


r = - --J'.P  p 
2 d Jin  n 


(5.4) 


it  represents  the  attenuation  of  the  front  pressure  in  the 
attenuation  region. 


5.2  Behaviour  of  Pn 

As  mentioned  in  the  last  section  PQ  is  the  front  overpressure 
at  the  beginning  of  the  attenuation  region  in  position  n = 16 . 
The  values  of  PQ , obtained  from  the  curve  fitting  based  on 
pressure  measurements  from  position  n - 16  to  the  tube  outlet, 
represent  fairly  well  the  measured  values  of  front  pressure 
at  n = 16. 

PQ  is  determined  by  the  events  in  the  build-up  region.  Since f 
however,  the  build-up  region  has  been  defined  to  be  the  16 
tube  diameters  closest  to  the  charge,  PQ  is  also  to  some 
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extent  determined  by  the  effects  which  are  characteristic 
for  the  propagation  of  the  shock  front  in  the  attenuation 
region.  This  is  for  example  seen  from  the  dependence  of 
PQ  on  wall  roughness.  Pq  decreases  with  increasing  roughness. 

We  will  more  easily  see  how  PQ  qualitatively  depends  on  q 
and  D if  we  determine  the  parameters  a and  b so  that  a curve 

Po  = aqb  (5.5) 

is  fitted  to  the  experimental  data. 

The  calculated  values  of  a and  b together  with  the  correlation 
coefficient  R are  listed  in  table  5.2a  for  CFG  1 and  table  5,2b 
for  CFG  2. 


We  first  notice  the  relatively  good  correlation  for  this 
curve  form,  indicating  that  it  describes  fairly  well  how  PQ. 
increases  with  increasing  q. 

It  is  seen  that  b < 1.  This  means  that  an  increase  in  the 
charge  group  does  not  result  in  a proportional  increase  of 
Pq  * 

We  also  notice  as  a trend  that  a increases  with  increasing  D, 
while  b decreases  with  increasing  D.  This  should  mean  that 
PQ  is  larger  for  the  larger  scales,  but  that  the  difference 
between  the  scales  is  smaller  for  the  larger  charge  groups. 
The  trends,  however,  are  weak,  and  one  should  be  careful  in 
drawing  conclusions. 


5.3  Attenuation  region 


In  this  section  the  attenuation  of  the  front  pressure  will 
be  discussed,  that  is,  the  dependence  of  the  attenuation  on 
scale,  charge  group  and  roughness.  Since  C2  represents  the 
attenuation  of  the  front  pressure,  we  shall  study  how  C?. 
depends  on  these  parameters.  5 
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We  shall  not  make  any  quantitative  statements,  but  concentrate 
on  some  qualitative  features.  We  shall  do  this  by  fitting 
simple  curve  forms  to  the  values  found  for  C2.  That  makes 
it  easier  to  see  qualitatively  how  C-2  depends  on  the  different 
parameters.  It  should  not  be  understood  as  an  attempt  to 
describe  the  real  behaviour  of  C2. 

5.3.1  Dependence  of  C2  on  scale 

We  shall  first  look  at  the  dependence  of  C2  on  tube  diameter  D. 
The  values  of  C2  and  D are  used  to  determine  a linear  curve 

c2  = a + bD  ( 5 . G ) 

and  the  results  are  listed  in  table  5.3a  for  CFG  1 and 
table  5.3b  for  CFG  2. 

In  CFG  1,  b < 0 in  all  cases,  except  for  one.  This  means 
that  C2  decreases  when  D increases,  showing  that  the  attenu- 
ation of  the  front  pressure  as  the  shock  propagates  through 
the  tube  is  stronger  for  the  smaller  scales. 

In  table  5.4a  a linear  curve  in  charge  group  q is  fitted  to 
the  values  of  b.  The  tendency  is  that  b decreases  when  q 
increases.  Since  b measures  the  difference  in  attenuation 
between  the  three  scales,  it  means  that  this  difference 
becomes  larger  for  larger  charge  groups. 

For  CFG  2 , b < 0 in  8 cases  and  b > 0 in  7 cases.  That  is, 
no  tendency  of  the  attenuation  to  decrease  with  increasing 
scale  is  seen.  This  difference  between  CFG  1 and  CFG  2 might 
be  explained  statistically,  but  it  may  also  result  from  the 
different  physical  situations  in  the  two  configurations , as 
will  be  discussed  in  section  5.7. 

From  table  5.4b,  where  a linear  curve  in  q is  fitted  to  b, 
we  see  that  for  sufficiently  large  charge  groups,  b should 
become  negative,  as  in  CFG  1. 
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5.3.2  Dependence  of  C2  on  charge  group 
The  curve 

c2  = a q°  (5.7) 

is  fitted  to  the  values  of  C2.  The  results  are  listed  in 
table  5.5a  for  CFG  1 and  table  5.5b  for  CFG  2.  In  both 
configurations  b is  negative,  so  that  C2  decreases  when  q 
increases.  It  means  that  the  attenuation  is  smaller  for  the 
larger  charge  groups . 

We  also  notice  a weak  tendency  of  b to  decrease  with  increasing 
D.  If  that  is  true,  it  should  mean  that  the  effect  of  smaller 
attenuation  for  larger  charge  groups,  is  increasing  with  in- 
creasing scale.  This  is  consistent  with  what  was  mentioned 
in  the  previous  section. 

It  is  a trend  that  b increases  with  increasing  value  of  rL, 
which  would  mean  that  the  presence  of  large  scale  wall  rough- 
ness tends  to  oppose  the  decreasing  effect  of  increasing  charge 
group  on  the  attenuation.  In  other  words,  the  damping  of  the 
front  pressure  when  large  scale  wall  roughness  is  present,  is 
not  so  much  decreased  by  increasing  charge  group,  as  it  other- 
wise would. 

<1 

5.3.3  Dependence  of  C2  on  large-scale  wall  roughness 

As  expected  C2  increases  with  increasing  value  of  rL.  The 
attenuation  of  the  front  pressure  is  larger  with  larger  rough- 
ness present. 

In  tables  5.6a  for  CFG  1 and  table  5„6b  for  CFG  2,  the  linear 
form 


C2  = a + b rp 

is  fitted  to  the  values  of  C2 . 


(5.8) 
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b is  positive,  so  that  C2  increases  with  increasing  r^,. 
b decreases  weakly  with  increasing  scale.  The  increase  in 
the  attenuation  with  increasing  rL  is,  in  other  words,  a 
little  stronger  for  the  smaller  scales.  The  damping  effect 
of  a large  wall  roughness  is  thus  somewhat  stronger  for  the 
smaller  scales  than  the  larger  ones. 

With  the  presence  of  a large  scale  roughness  r^  >0,  a deviation 
from  scaling  with  respect  to  dimension  is  introduced. 


One  should,  however,  be  a little  careful  in  drawing  conclusions. 
The  decrease  in  b with  scale  is  weak,  and  occurs  in  neither 
of  the  configurations  between  all  the  three  scales  for  all 
charge  groups . 

5.3.4  Summary  of  the  curve  fitting  procedure 

The  discussion  in  the  preceeding  sections  can  be  summarized 
by  fitting  the  linear  expression 


C2  = ^1  + ^2  D + 63  rL  + ^4  *3 

+ 3s  D rL  + 3g  D £n  q + g_  rL  in  q 


(5.9) 


to  the  values  of  C2  for  both  configurations.  Notice  that  the 
dependence  of  C2  on  charge  group  q in  this  expression  is 
through  in  q. 

The  resulting  values  of  6 and  the  partial  correlations  are 
listed  in  table  5.7.  Some  of  the  partial  correlation  coeffi- 
sier.ts  are  not  significant.  The  scatter  of  the  values  of  C2 
is  too  large  to  expect  a good  correlation  to  eq  (5.9). 

However,  eq  (5.9)  is  not  expected  to  give  a quantitative 
description  of  the  dependence  of  C2  on  q,  rL  and  D,  but  only 
to  reveal  some  of  the  qualitative  features  of  this  dependence. 
To  decide  if  the  dependence  of  C2  on  each  of  the  terms  in 
eq  (5.9)  is  significant,  one  should  not  only  look  at  the 
partial  correlation  coefficients  in  table  5.7,  but  also  at 
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the  curve  fitting  procedures  in  the  preceeding  sections. 

With  the  values  of  (3-^  inserted  in  eg  (5.10)  , it  takes  the 
form 


C2  = 1.09  + 2.0  D + 27.6  rL  - 0.094  In  q 

- 20  D r^  - 0.55  D £n  q - 2.2  r^  £n  q 


(5.11a) 


for  CFG  1 and 


C2  = 1.47  + 2.75  D ~ 6.16  rL  - 0.175  £n  q 
- 24  D rjj  ~ 0.42  D £n  q + 2.7  rL  £n  q 


(5.11b) 


for  CFG  2. 

a2c2 

We  notice  that  -r-rr — — and 

3D  3rL 


32C< 


3D  3 £n 


q 


are  negative  in  both  configurations,  indicating  that  the 
increase  in  C2  with  increasing  scale  is  larger  for  larger 
charge  groups  and  larger  values  of  wall  roughness . 

32C2 

3 5 — p — has  different  sign  in  the  two  configuration. 


5 . 4 Comparisison  with  exponential  curve  form 

In  this  chapter  the  front  pressures  in  the  attenuation  region 
have  been  fitted  to  the  curve  form  eq  (5.1).  It  has  been 
proposed  that  the  decay  in  front  pressure  follows  an  expo- 
nential  law  (7),  that  is,  that  the  front  pressure  in  position 
n can  be  expressed  as 

-C2 (n-16) 

P = Po  e (5.12) 


where  PG  is  the  front  pressure  in  position  n = 16 . 
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We  have  not  used  this  function  for  curve  fitting  to  all 
our  data,  only  for  the  case  of  rL  = 0 in  CFG  1 we  have  made 
a comparison  between  the  two  curve  forms . The  comparison 
is  made  for  each  charge  group  and  scale  separately.  The 
method  used  is  one  described  by  T J Dylewski  in  ref  48)  . 

The  curVe  forms  are  fitted  to  the  data  and  the  square  devi- 
ation  S is  calculated  for  both  curves.  Then  all  the  pressures 
in  positions  n = 16  to  n = 32  are  fitted  to  both  curve  forms 
and  the  square  deviation  S'  is  calculated  for  both.  In  a 
similar  manner  the  square  deviation  S"-for  the  fitting  of 
the  pressures  in  positions  n = 36  to  n = 48  is  calculated 
for  both.  The  curve  form  giving  the  smallest  value  of 
g 

q-j — — g „ is  the  best  one. 

This  has  been  done  for  each  charge  group  and  scale  for  CFG  1, 
rL  = 0,  and  the  values  are  listed  in  table  5.8.  The  test 
gives  different  results  for  the  different  scales  and  charge 
groups , with  each  curve  form  preferred  about  the  equal 
numbers  of  times. 


5 . 5 Shots  with  small  charges 

The  front  pressures  measured  in  the  shots  with  the  small 
charges  in  scale  3:100  were  also  used  to  determine  Ci  and 
C2  in  eq  (5.1).  The  resulting  values  of  C± , C2,  PQ  and 
the  correlation  coefficient  R are  listed  in  tables  5.9a 
and  5.9b  for  CFG  1 and  CFG  2 respectively.  We  notice  the 
good  correlation  obtained  also  in  this  case. 

5.5.1  Build-up  region 

Tne  values  of  PQ  are  fitted  to  eq  (5.5)  giving 


P0  = 4.0  q 


0.64 


with  a correlation  R « 1.00  for  CFG  1,  and 
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PD  = 5.5  q 


0.54 


with  R = 0.93  for  CFG  2.  When  this  is  compared  with  the 
corresponding  results  for  the  larger  charges  in  the  main 
experimental  series  (see  tables  5.2a  & b) , 


P„  = 36.8  q 


0.15 


for  CFG  1 and 


P 


o 


28.2 


q0'20 


for  CFG  2,  it  is  seen  that  PQ  increases  faster  with  charge 
for  the  smaller  charges  than  for  the  larger.  The  increase 
in  PQ  is  not  proportional  to  the  increase  in  q,  however. 

PQ  as  a function  of  q is  shown  in  figure  '5.1a  for  CFG  1 
and  in  figure  5.1b  for  CFG  2.  The  figures  indicate  that 
there  is  a limit  on  the  pressures  which  can  be  obtained  by 
increasing  the  charge  weight. 


5.5.2  Attenuation  region 

A linear  expression  in  in  q is  fitted  to  the  values  of  C2 / 
resulting  in 


C2  = 1.56  - 0.198  £n  q 
with  correlation  R = 0.97  for  CFG  1 and 
C2  = 0.88  - 0.009  in  q 

with  the  poor  correlation  R = 0.14  for  CFG  2 

When  this  is  compared  with  the  corresponding  results  for 
the  main  series  (eg  5.11). 

C2  = 1.39  - 0.177  in  q 
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for  CFG  1 and 

C2  = 1.88  - 0.238  Jin  q 

for  CFG  2,  it  is  seen  that  the  decrease  in  C2  with  increasing 
charge  is  about  the  same  or  somewhat  larger  for  the  larger 
charges  than  the  smaller  ones. 

C2  is  shown  as  a function  of  q in  figure  5.2a  for  CFG  1 
and  figure  5.2b  for  CFG  2. 


5 . 6 Comparison  with  NDCS 

The  results  in  this  chapter  will  now  be  compared  with  results 
obtained  at  the  Norwegian  Defence  Construction  Service  (NDCS) 

(6) . Two  of  the  tubes  described  in  ref  (6)  were  of  steel. 

They  were  open  in  both  ends , and  the  charges  were  detonated 
in  the  middle  of  the  tubes.  The  tube  diameters  were  D = 0.05  m 
and  D = 0.20  m,  corresponding  to  scales  1:100  and  4:100, 
respectively.  The  large  scale  roughness,  as  defined  in  this 
report,  was  rL  = 0 . In  the  tube  with  D = 0.05  m the  charges 
used  ranged  from  9.5  g to  17 . 5 g , corresponding  to  charge 
groups  q = 76  kg/m3  and  140  kg/m3.  In  the  tube  with  D = 0.20  ra, 
the  charge  weights  were  in  the  range  9.5  g to  200  g,  corre- 
sponding to  charge  groups  q = 1.1875  kg/m3  to  25  kg/m3. 

Sq  (5.1)  was  fitted  to  the  front  pressures,  resulting  in 

C-,  = 18 

JL 


and 


C2  = 0.66 

These  values  are  in  good  accordance  with  our  results. 

C2  - 0.66  is  larger  than  the  values  of  C2  in  our  main  series 
of  shots,  but  are  in  the  range  of  values  of  C2  in  the  shots 
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with  small  charges.  The  q-values  in  our  shots  with  small 
charges  are  included  in  the  range  of  q-values  in  the  NDCS 
experiments.  Also  Cg  = 18  is  smaller  than  the  mean  of  the 
Cg-values  in  our  main  series,  and  is  closer  to  the  Cg-values 
obtained  from  the  shots  with  small  charges. 


5.7  Difference  between  CFG  1 and  CFG  2 


In  CFG  1 the  tube  end  closest  to  the  charge  is  open  while 
in  CFG  2 it  is  closed.  This  gives  rise  to  important  diffe- 
rences in  the  blast. 

In  CFG  1 a rarefaction  wave  will  be  reflected  from  the  open 
end  and  propagate  in  the  heated  gas  behind  the  shock  front 
with  a larger  speed  in  the  same  direction.  Sooner  or  later 
it  will  overtake  the  shock  and  attenuate  it.  The  same  will 
happen  in  CFG  2,  but  the  wave  reflected  from  the  closed  end 
will  be  a shock  wave,  thus  increasing  the  front  pressure  of 
the  first  shock  when  it  sooner  or  later  catches  it  up. 

As  a rule  of  thumb  the  front  pressure  in  CFG  2 should  be 
the  same  as  the  front  pressure  in  CFG  1 resulting  from  a 
charge  twice  as  large. 

* 

What  we  observe,  however,  is  smaller  front  pressures  in 
CFG  2 than  in  CFG  1 in  the  tube.  This  must  result  from  the 
other  difference  between  the  two  configurations.  In  CFG  1 
the  charges  are  detonated  in  the  tube,  while  in  CFG  2 they 
are  detonated  in  a detonation  chamber  with  a cross-sectional 
area  twice  as  large  that  of  the  tube.  When  the  energy  from 
a detonation  in  CFG  2 propagates  from  the  chamber  into  the 
tube,  it  must  be  delayed  by  the  decrease  in  cross-sectional 
area. 

One  difference  between  CFG  1 and  CFG  2 can  be  seen  from  the 
curve  fitting  in  section  5.2.  Tables  5.2a  & b give  the 
values  of  a and  b in  the  equation. 
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t>  b 

po  = a q (5.5) 

It  is  seen  that  b is  larger  for  CFG  2 than  for  CFG  1 , and 
that  a is  smaller  for  CFG  2 than  for  CFG  1.  Thus,  as  seen 
from  Tables  5.2,  the  values  of  PQ  for  q = 80  kg/m3  are 
larger  for  CFG  1 than  for  CFG  2,  but  the  increase  with 
increasing  q is  greater  in  CFG  2. 

We  have  not  studied  systematically  the  difference  in  pressure 
time  histories  of  the  blast  in  the  two  configurations,  but 
a short  examiation  seems  to  indicate  a difference  as  shown 
in  figure  5.3.  The  blast  has  a larger  front  pressure,  but 
a shorter  positive  duration  in  CFG  1 than  in  CFG  2.  In  CFG  2 
a second  peak  is  following  the  first  peak  some  distance 
behind  without  overtaking  it  during  the  propagation  through 
the  tube . 


6 DISCUSSION 

6 . 1 Build-up  region 

In  section  5.2  it  was  seen  that  PG , the  front  pressure  in 
position  n = 16  at  the  beginning  of  the  attenuation  region, 
is  an  increasing  function  of  charge  group,  but  that  the 
increase  was  smaller  than  a proportional  increase.  It  was 
also  seen  that  PQ  increases  with  increasing  scale,  but  that 
this  increase  is  smaller  for  the  larger  charge  groups.  To 
that  extent  PG  obtained  from  the  curve  fitting  really  is 
determined  by  the  build-up  phase  of  the  shock,  it  seems  to 
indicate  that  there  is  an  upper  limit  on  the  front  pressures 
which  can  be  obtained  by  increasing  the  loading  density. 

This  would  be  in  accordance  with  a discussion  by  D R Curran 
in  ref  (9) . He  states  that  a wide  range  of  initial  pressure 
in  the  detonation  chamber  will  produce  initial  air  shock 
pressures  which  lie  close  together  in  a region  about  100  bars. 
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We  shall  not  go  further  into  this  problem.  It  would  have 
been  or  interest,  however,  to  study  not  only  the  front  pres- 
sure, but  the  whole  pressure-time  history  of  the  blast  at 
the  measurement  stations  closest  to  the  charge. 

Since  it  seems  to  be  an  upper  limit  on  the  initial  front 
pressures,  it  must  also  be  an  upper  limit  on  the  deviation 
from  scaling  with  respect  to  dimension.  That  is,  the  initial 
front  pressure  can  not  be  very  much  larger  in  a full  scale 
test  than  in  a small  scale  test.  So  for  the  problem  of 
extrapolating  from  model  tests  to  full  scale  tests  it  seems 
to  be  more  important  to  concentrate  on  the  attenuation  of  the 
front  pressure  of  the  shock  than  its  initial  formation. 


6. 2 Attenuation  region 


There  are  two  main  sources  for  the  attenuation  of  the  front 
pressure  as  the  shock  propagates  through  the  tube  (6) , (i) 

rarefaction  waves  which  gradually  overtake  and  eat  up  the 
shock  front,  (ii)  wall  roughness. 

We  defined  a large-scale  wall  roughness  r-^  in  section  3.1, 
eq  (3.1)  and  our  tubes  provide  three  degrees  of  large  scale 
roughness,  r-^  = 0,  r^  = 0.03  and  rp  = 0.06.  This  roughness 
leads  to  a complicated  reflection  and  diffraction  pattern 
of  the  shocks  (10) , giving  rise  to  an  attenuation  of  the 
fronu  pressure. 

There  is  an  additional  roughness  which  we  shall  call  the 
small  scale  roughness,  rg.  The  walls  in  the  tubes  are  not 
quite  smooth.  Both  with  rL  = 0,  rL  = 0.03  and  rL  = 0.06 
there  are  irregularities  of  the  size  of  grains  of  sand  in 
tne  steel  walls.  We  have  measured  the  mean  height,  k,  of 
these  irregularities.  For  all  the  tubes  in  scale  1:100  it 
was  found  to  be 


k = 0.044  mm 


(6.1a) 
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and  for  all  the  tubes  in  scales  2:100  and  3:100  it  was 


k = 0.080  mm 


(6.1b) 


We  shall  define  the  small-scale  roughness  rg  by  the  relation 
k 


rS  “ d 


(6.2) 


This  gives 
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The  total  pressure  loss  can  be  expressed  by  the  relation 


dP 

dn 


(|El  + 

rdP' 

+ 

fdPl 

Mr 

dn_ 

rL 

dn 

s / 

(6.4) 


where 


dP 


dn 


R 


is  the  pressure  loss  due  to  overtaking  rare- 


faction waves , and 


dP' 

and 

dP' 

dn_ 

L 

rs 


are  the  pressure  losses 


due  to  large-scale  and  small-scale  roughness,  respectively. 


An  understanding  .of  each  of  these  damping  effects  and  of 
their  relative  importance  would  give  an  understanding  of 
the  total  pressure  attenuation,  and  then  also  of  the  depen- 
dence of  the  attenuation  on  scale. 


6.3  Pressure  loss  due  to  rarefaction  waves 


Many  models  have  been  proposed  for  the  shockfront  attenuation 
due  to  overtaking  rarefaction  waves.  One  of  them,  proposed 
by  Porzel  (11)  states  that 


f (P) 


1 3P 
U 3 1 


(6.5) 
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U is  the  shock  velocity,  is  the  slope  of  the  rarefaction 

wave  and  f (P)  is  a function  of  the  front  overpressure. 

In  an  air  driven  shock  tube,  where  ~ = 0,  there  is  no 

0 t, 

pressure  loss  due  to  rarefaction  waves.  For  a peaked  blast 
3P 

profile,  is  large,  and  the  pressure  loss  will  be  large. 

We  have  not  examined  the  pressure-time  profiles  from  our 
experiments , and  consequently  we  can  not  determine  the 
relative  importance  of  the  pressure  loss  due  to  rarefaction 
waves.  Also,  a detailed  study  of  the  difference  in  pressure- 
time profiles  between  CFG  1 and  CFG  2 would  give  information 
about  the  rarefaction  waves,  since  the  tube  end  close  to  the 
charge  is  open  in  CFG  1 and  closed  in  CFG  2. 


6 . 4 Pressure  loss  due  to  large-scale  roughness 

We  found  in  section  5.3.3  that  the  increase  in  C2  with 
increasing  rL  is  larger  for  the  smaller  scales,  which  means 
that  the  damping  effect  of  a non  zero  rj,  is  larger  in  the 
smaller  scales.  A scale  dependence  of  the  front  pressure 
therefore  results  from  the  presence  of  a large-scale  rough- 
ness. 

* 

The  value  of  C2  is  high  also  in  the  case  of  rL  = 0,  indicating 
that  the  large  scale  roughness  can  not  be  the  only  important 
source  of  the  attenuation.  A scale  dependence  of  Cp  is  also 
present  for  rL  = 0 , so  there  must  be  additional  effects  which 
do  not  scale. 

6 . 5 Pressure  loss  due  to  small-scale  roughnes s 

From  eq  (6.3)  it  was  seen  that  the  small  scale  roughness  r<- 
in  our  experiments  does  not  scale f i e,  it  has  not  the  same 
values  in  all  the  three  scales*  To  the  extent  the  pressure 
loss  due  to  this  roughness  is  important  for  the  total  pressure 


\ 
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loss,  it  will  lead  to  a dependence  of  the  attenuation  on 
scale . 

The  values  of  C2  are  relative  large  even  in  the  case  of 
rL  = 0 , indicating  that  if  the  assumption  in  eq  (6.4)  is 
valid,  an  important  part  of  the  attenuation  must  be  due  to 
the  small  scale  roughness  and/or  rarefaction  waves.  Also, 
since  C2  is  scale  dependent  even  when  rL  = 0 , there  must,  be 
a scale  dependence  in  at  least  one  of  these  effects.  Since 
we  have  not  studied  the  effect  of  the  rarefaction  waves,  we 
can  not  decide  if  this  gives  rise  to  a scale  dependence  or 
not.  Because  of  the  difference  between  CFG  1 with  the  open 
tube  end  and  CFG  2 with  the  closed  end  with  respect  to  rare- 
faction  waves , one  should  perhaps  expect  a larger  difference 
in  the  values  of  C2  in  the  two  configurations  than  we  have 
observed,  if  overtaking  rarefaction  waves  should  be  the  main 
reason  for  the  front  pressure  attenuation  when  rL  = 0 . 

Several  experimental  and  theoretical  investigations  of  the 
front  pressure  attenuation  in  air  driven  shock  tubes  indicate 
that  a small  scale  roughness  can  result  in  an  appreciable 
front  pressure  attenuation  as  a shock  propagates  through 
the  tube  (12),  (13),  (14),  (15).  The  magnitude  of  the  attenu- 

ation over  a distance  50  to  100  tube  diameters  may  be  in  the 
range  of  20  to  80  percent. 

t 

We  shall  study  the  possibility  that  the  main  part  of  the 
deviation  from  scaling  with  respect  to  dimension  when 
rL  = 0 might  be  explained  by  the  pressure  loss  due  to  small- 
scale  roughness. 


6 . 6 The  physical  reason  for  the  pressure  loss  due  to  small- 
scale  roughness 

The  pressure  loss  due  to  the  small-scale  roughness  may  be 
explained  in  the  same  way  as  for  the  large  scale  roughness, 
namely  that  the  projections  give  rise  to  a set  of  reflections 
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and  diffractions  of  the  shock  which  attenuates  the  front 
pressure  (10) . It  may  also,  however,  result  from  the  viscous 
friction  with  the  tube  wails,  arising  in  a turbulent  boundary 
layer  behind  the  shock. 

It  is  difficult  to  decide  which  of  these  explanations  is 
the  most  important  one  without  a quantitative  and  more  de- 
tailed study  of  the  attenuation  observed  in  our  experiments. 

We  shall  in  the  following  present  a possible  formulation  of 
the  problem  if  the  pressure  loss  is  caused  by  viscous  friction. 


6.7  Pressure  loss  due  to  viscous  friction 


The  pressure  loss  dP  due  to  viscous  friction  over  a tube 
section  of  length  dn  is  assumed  to  be  independent  of  position 
n and  proportional  to  the  pressure  at  this  position,  i e: 


dP 

dn 


KP 


where  K is  a constant. 


(6.6) 


Ey  the  Rankine-Hugoniot  equations,  this  can  also  be  expressed 
as 


^ = K (1  + y2)  (MS2  - 1) 


(6.7) 


where  Ms  is  the  Mack  number  of  the  shock  front,  or,  for 
strong  shocks,  approximately  as 


dP 

dn 


(1  + y2 

1 " M2) 


(6.8) 


where  Mf  is  the  Mach  number  of  the  flow  behind  the  shock. 
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Eq  (6.6)  was  proposed  by  Emrich  and  Curtiss  (16) 
after  experiments  with  air  driven  shock  tubes  of  different 
size  and  shape.  Their  tubes  had  zero  large-scale  roughness, 
and  no  rarefaction  waves  were  present.  They  found  K to  be 
almost  independent  of  n and  only  slowly  varying  with  P. 

The  value  of  K is  about  3.5  • 10  Their  experiments 

comprised,  however,  only  shock  strengths  up  to  Mg  = 3. 

dP 

One  must  expect  ^ to  depend  on  the  whole  pressure-time 
history  of  the  blast  when  passing  position  n.  Specially, 
it  must  depend  on  the  positive  duration  t+  of  the  blast. 

At  BRL  K was-  proposed  to  be  of  the  form  (3) 

I<  = K'  (a1  + 2_)  (6.9) 

with  K'  and  a'  approximately  constant. 

Eq  (6.8)  is  of  the  same  form  as  the  standard  formula  for 
the  pressure  loss  in  a subsonic  flow  in  a tube, 

dP  r i 2i 

dH  = M*  P w J (6.10) 

where  w is  the  flow  velocity  and  p the  fluid  density. 

X is  a factor  which  depends  on  the  small  scale  roughness 
hnd  the  Reynolds  number.  For  large  Reynolds  numbers  X 
becomes  independent  of  Reynolds  numbers,  and  can  be  expressed 
by  the  Nikuradse  formula 


. _ 0.25 

(l.g  3/715x2  (6.11) 

rs 

lg  denotes  the  logarithm  with  base  10. 

We  shall  assume  that  the  small-scale  roughness  enters  into 
eq  (6.6)  for  a blast  wave  through  the  same  factor  as  in 
eq  (6.10)  for  a subsonic  flow,  and  since  the  Reynolds  numbers 
-n  our  case  are  sufficiently  large, the  Nikuradse  formula 
eq  (6.11)  will  be  used. 
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We  shall  in  other  words  assume  that  K in  eqs  (6.6  - 6.8) 
is  of  the  form 


K 


a 


0.25 


tig  (—)) 


= aX 


(6.12) 


rs 


where  a according  to  eq  (6.9)  depends  on  the  positive  duration 
t+  of  the  blast. 


6 • 8 Magnitude  of  pressure  loss  due  to  viscous  friction 

If  we  knew, the  formulae  describing  the  pressure  losses  in 
eq(6.4)  this  equation  could  be  integrated  to  give  the  front 
pressure  as  a function  of  position  in  the  tube  for  the  various 
tube  dimensions,  charge  weights  and  values  of  large-scale 
roughness.  However,  this  study  will  be  restricted  to  a sug- 
gested relation  for  the  pressure  loss  due  to  small-scale 
roughness  when  it  is  assumed  to  result  from  viscous  friction 
to  the  walls. 

We  shall  examine  if  this  deviation  due  to  small-scale  rough- 
ness can  be  responsible  for  the  total  deviation  from  scaling 
in  the  attenuation  of  the  front  pressure.  That  is,  we  must 
examine  if  there  is  a value  of  a in  eq  (6.12)  which  is  suf- 
ficiently large  to  be  responsible  for  all  the  deviation  ob- 
served and  sufficiently  small  to  avoid  an  unrealistic  large 
pressure  loss  due  to  the  small-scale  roughness.  We  shall  do 
this  by  trying  to  determine  the  values  P'  which  the  front 
pressure  would  have  had  at  different  positions  in  the  tube 
for  the  three  scales,  if  there  had  been  no  pressure  loss  due 
to  snail  scale  roughness.  If  all  the  deviation  should  be 
explained  by  this  effect,  the  attenuation  of  P'  should  be 
the  same  for  all  scales. 

Since  this  is  more  an  illustration  than  a systematic  analysis, 
we  shall  do  it  only  for  r^  -•  0 in  CFG  1. 
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The  pressure  loss  AP  over  a tube  length  of  4 tube  diameters 
between  two  succeeding  measurement  stations  in  positions, 
n-4  and  n are  calculated  according  to  eqs  (6.7)  and  (6.12). 


AP  = 4‘|*aX (Ms2  - 1)  (6.13) 

B D 

The  value  of  the  shock  velocity  U is  set  equal  to  where 
At  is  the  difference  in  arrival  times  between  the  positions 
n+4  and  n-4.  The  value  of  U obtained  in  this  way  is  asso- 
ciated with  position  n.  Values  of  U are  shown  in  table  6.1 
for  rx,  = 0,  CFG  1. 


P'  is  obtained  at  a measurement  station  in  position  n by 

adding  to  the  measured  pressure  P the  sum  of  all  the  pressure 

losses  between  previous  measurement  stations.  Because  of 

. 

the  great  scatter  of  P,  we  shall  use  the  values  P = C]_ 
from  the  curve  fitting  in  chapter  5. 
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The  arguments  presented  are  of  a somewhat  approximate  chara- 
cter, and  a value  of  a can  not  be  determined  exactly.  It  is 
determined  by  the  condition  that  the  attenuation  of  P!  should 
be  independent  of  scale.  If  the  attenuation  is  independent 


of  scale,  then  the  values  of 


P' (3:100) 


and 


P' (2:100) 


are 


P' (2:100)  P' (1:100) 

invariant  with  regard  to  n,  and  they  are  equal  to 

Pn  (3:100)  , Po(2:100)  . . . , _ . . 

P T -TOO)  and  P (!"• Too)  respectively,  where  P0  denotes  an 

initial  pressure.  This  turns  out  to  be  approximately  satis- 
fied for  a = 0.6  as  is  seen  in  table  6.2. 


Figure  6.1  shows  an  example  of  P'  as  function  of  n together 
with  the  values  of  P obtained  from  the  curve  fitting  for  the 
three  scales. 

In  position  n = 48  the  pressure  loss  due  to  eqs  (6.6)  and 
(6.12)  with  a = 0.6  is  about  60%  of  the  total  pressure  loss.  It 
then  follows  that  if  the  pressure  loss  from  the  small  scale 
roughness  can  be  responsible  for  about  60%  of  the  total 
pressure  loss  as  the  shock  propagates  through  the  tube,  it 
can  also  be  responsible  for  the  deviation  from  scaling  ob- 
served in  the  attenuation  of  the  front. 


60%  is  a very  likely  value  for  the  small  scale  roughness 
pressure  loss.  Based  on  experiments  at  BRL  a value  between 
20%  and  80%  has  been  proposed  (13) . Bobin  and  Thery  (12) 
have  done  an  examination  of  the  pressure  attenuation  with 
help  of  a computational  program.  Although  the  pressures  they 
have  examined  are  much  smaller  than  ours,  and  the  tubes  are 
much  longer,  their  results  should  indicate  a value  of  about  50%. 

As  a conclusion,  it  seems  to  be  a possible  hypothesis  that 
the  scale  dependence  in  front  pressure  attenuation  when 
rL  = 0,  may  result  from  viscous  friction  due  to  small-scale 
roughness . 

It  is  emphasized , however,  that  to  prove  this  hypothesis 
it  is  necessary  to  study  the  experimental  data  in  much 
further  detail. 


Summary 

Small-scale  roughness,  large-scale  roughness  and  rarefaction 
waves  are  supposed  to  be  the  main  reasons  for  the  attenuation 
of  the  front  pressure.  The  data  have  not  been  sufficiently 
analysed  to  give  any  quantitative  statements  of  the  relative 
importance  of  each  of  these  effects.  It  is  reasonable  to 
assume  that  the  deviation  from  scaling  with  respect  to  dimen- 
sion is  mainly  due  to  the  small  scale  roughness.  It  is  not 
understood  if  the  small-scale  roughness  pressure  loss  results 
primarily  from  viscous  friction  or  from  the  diffraction  pattern 
of  the  shock  from  the  small  wall-irregularities.  Our  data 
have  not  been  analysed  in  such  a way  that  we  can  decide  which 
of  these  are  most  important.  If  the  damping  is  assumed  to 
result  mainly  from  viscous  friction,  it  is  sufficient  that 
this  damping  is  about  60%  of  the  total  damping  at  the  end 
of  the  tube,  to  associate  all  the  scale  dependence  of  the 
attenuation  with  the  damping  from  small-scale  roughness  rs. 
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7 EXTRAPOLATION  TO  FULL  SCALE 

7 . 1 Effect  of  small-scale  roughness 

If  the  assumptions  in  chapter  6 are  valid,  the  front  pressures 
in  a full  scale  steel  tube  with  D = 5 m and  small  scale  rough- 
ness rg  can  easily  be  found  by  using  the  values  of  P'  obtained 
from  model  tests  and  subtracting  the  pressure  loss  due  to 
small  scale  roughness  as  given  by  eq  (6.6), 

An  example  of  full  scale  pressures  obtained  this  way  in  steel 

-5 

tube  with  k = 0.080  mm  (rg  = 1.6  • 10  ) is  shown  as  the  full 
drawn  line  in  figure  6.1.  Thismeans  that  if  the  assumptions 
in  chapter  6 are  valid  the  shock  parameters  in  a full  scale 
steel  tube  can  be  predicted  by  experiments  in  small  scale 
steel  tubes. 

Figure  7.1  shows  as  an  example  how  the  front  pressure  in 
position  n = 48  at  the  end  of  the  attenuation  region  depends 
or.  scale,  according  to  the  formulae  in  chapter  6.  k = 0.080  mm 
is  used. 

7 . 2 Real  ammunition  storage  site 

4 

Real  undergound  ammunition  storage  sites  are  not  made  of 
steel,  but  are  excavated  in  rock.  For  predicting  shock  para- 
meters in  this  case  from  the  shock  parameters  in  the  full 
scale  steel  case  it  is  necessary  with  a knowledge  of  the 
different  influence  of  steel  and  rock  on  the  blast.  Such 
an  analysis  will  not  be  attempted  here. 

It  is  also  necessary  to  determine  the  mechanism  which  causes 
the  attenuation  of  the  front  pressure  due  to  the  actual  rough- 
ness in  the  tunnel  walls.  The  front  pressures  obtained  in 
a full  scale  steel  tube  must,  however , be  assumed  to  represent 
an  upper  bound  for  the  pressures  which  can  be  obtained  in  a 
geometrically  similar  tube  in  rock. 
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The  geometrical  conditions  in  an  actual  storage  site  are 
not  as  simple  as  in  our  models.  We  have  not  examined  the 
effects  on  the  blast  resulting  from  more  complex  configura- 
tions. Effects  not  encountered  in  our  simple  models  may 
enter  into  the  problem  under  more  complicated  geometrical 
conditions . 


7.3  Full  scale  tests 


It  would  be  of  great  importance  to  compare  data  from  mcdel- 
and  full  scale  experimental  series  from  several  institutions 
to  single  out  how  different  effects  scale.  We  have  not  done 
such  an  examination  of  all  existing  experimental  information. 
We  have  only  studied  the  results  from  some  full  scale  tests 
performed  in  Norway,  and  compared  them  with  the  present  model 
tests.  There  are  very  many  differences  between  the  models 
and  these  full  scale  tests,  so  there  are  severe  limitations 
on  the  conclusions  which  can  be  drawn. 

A short  review  of  the  full  scale  model  tests  will  be  given 
in  the  rest  of  this  chapter. 

7.3.1  Tests  at  Raufoss 

The  Raufoss  experiments  are  described  in  ref  (17) . The 
tunnel  was  closed  in  one  end,  and  the  diameter  was  2.65  m, 
corresponding  to  scale  53:100.  In  some  shots  the  charges 
were  detonated  in  a sidechamber,  while  in  four  shots  the 
charges  were  detonated  in  the  tunnel.  We  shall  look  at  the 
results  from  this  last  case.  They  must  be  compared  with  the 
results  from  shots  with  charges  twice  as  large  as  in  our 
CFG  I.  In  two  shots  the  charge  weight  was  1000  kg,  the  only 
difference  being  that  the  charge  in  one  case  was  closer  to 
the  tunnel  end  than  in  the  other.  The  charge  weights  in  the 
two  other  shots  were  300  kg  and  100  kg.  This  corresponds 
to  charge  groups  53.7  kg/m“ , 16.1  kg/m  and  5.37  kg/m  . 

5.37  kg/m  is  much  lower  than  we  have  had  while  the  other 
two  are  in  the  range  of  values  of  our  shots.  The  pressures 
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obtained  in  the  Raufoss  experiments  are  listed  in  table  7.1 
together  with  some  actual  pressure  measurements  from  our 
model  tests.  In  the  shots  with  small  charges  only  tubes  with 
rp  = 0 were  used. 

We  see  that  the  Raufoss  pressures  are  very  low.  This  may 
be  due  to  the  tunnel  roughness  and  to  the  fact  that  the  tunnel 
walls  are  of  rock  and  not  of  steel.  It  may  also  be  due  to 
the  considerable  geometrical  differences. 

We  have  fitted  eq  (5.1)  to  the  Raufoss  pressures  to  see  how 
the  attenuation  of  the  front  pressure,  represented  by  C2 , 
in  this  full  scale  case  deviates  from  the  attenuation  in  the 
model  tests. 

In  section  5.3  the  values  of  C2  from  CFG  1 in  the  main  series 
of  the  model  tests  were  represented  by  the  equation 

C2  = 1.09  + 2.0  D + 27.6  rL  - 0.094  £n  q 

(5.11a) 

- 20  D rp  - 0.55  D £n  q - 2.2  rp  In  q 

If  we  insert  in  this  equation  rL  = 0.06,  which  should  be 
the  most  probable  value  of  the  wall  roughness  in  the  Raufoss 
tunnel,  and  for  q two  times  the  actual  q-values  in  the  Raufoss 
shots,  we  can  compare  the  attenuation  found  in  the  Raufoss 
shots  with  the  attenuation  found  in  the  model  tests.  The 
results  are  shown  in  table  7.2.  They  seem  to  support  the 
conclusion  based  on  the  model  tests , that  the  attenuation 
of  the  front  pressure  decreaces  with  increasing  scale.  For 
the  two  Raufoss  shots  with  q = 53.7  kg/m  the  values  of  C2 
are  about  as  large  or  somewhat  larger  than  the  corresponding 
values  of  C2  from  the  model  tests.  This  may,  however,  be  due 
to  the  different  attenuating  effect  of  rock  and  steel. 


45 


7.3.2  Tests  at  Lista 

The  Lista  experiments  are  described  in  ref  (9) . The  tunnel 
was  open  in  both  ends  and  had  a square  cross  section.  The 
area  of  the  cross  section  was  about  12  m3 , which  corresponds 
to  a hydraulic  diameter  D = 4 m or  scale  30:100.  The  charges 
were  detonated  in  the  tunnel,  so  it  corresponds  to  our  CFG  1. 
The  shots  were  fired  with  four  charge  weights,  10  kg,  57  kg, 

166  kg  and  500  kg,  corresponding  to  charge  groups  0.156  kg/m3, 
0.89_  kg/m3,  2.59  kg/m3  and  7.81  kg/m3,  respectively . 

The  charge  groups  are  far  below  the  range  of  our  experiments. 
Consequently,  no  direct  comparison  of  the  measured  pressures, 
shown  in  table  6.3,  and  the  pressures  from  the  model  tests 
can  be  done. 

Eq  (5.1)  is  fitted  also  to  the  Lista-pressures  to  compare 
the  values  of  C2  such  obtained  with  the  values  of  C2  from 
-he  model  tests.  When  rL  = 0.06  as  the  most  probable  value 
of  wall  roughness  in  the  Lista  tunnel,  and  the  actual  q-values 
from  the  Lista  shots,  are  inserted  in  eq  (5.11a)  we  can  compare 
-he  attenuation  in  the  Lista  experiments  with  the  attenuation 
found  in  the  model  tests. 

The  values  of  C2  calculated  from  eq  (5.11a)  and  from  the  Lista 
shots  are  shown  in  table  7.4.  It  should  be  noticed,  however, 
that  to  use  eq  (5.11a)  with  these  q-values  is  a considerable 
extrapolation  from  the  q-values  used  to  determine  the  para- 
meters in  the  equation. 

In  spite  of  this,  also  the  Lista  experiments  seem  to  support 
the  conclusion  that  the  attenuation  of  the  front  pressure 
decreases  with  increasing  scale. 

7.3.3  Summary 

A comparison  between  model  tests  and  geometrically  and  other- 
wise different  full  scale  tests  seems  to  give  most  informa- 
tion when  the  attenuation  of  the  front  pressure  is  examined. 
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^2  seems  to  decrease  also  when  the  scale  increases  towards 
full  scale.  This  is  in  accordance  with  the  observations  done 
in  the  small  range  of  scales  1:100  - 3:100. 

8 CONCLUSION 

The  primary  objective  of  this  work  has  been  to  examine  the 
validity  of  scaling  laws  (eg  2.7  and  2.11)  for  explosions 
-n  underground  installations.  The  results  have  shown  that 
within  the  range  of  our  small  scales  the  deviation  from 
scaling  is  highly  significant.  Further , it  is  clear  that 
if  the  deviation  from  scaling  with  respect  to  dimension  should 
continue  to  increase  to  the  same  degree  outside  our  scaling 
range,  then  small  scale  testing  would  be  quite  unsuitable 
for  predicting  explosion  effects  in  the  full  scale  case. 

This  leads  up  to  the  secondary  objective,  to  find  possible 
corrections  to  the  scaling  laws , which  make  them  suitable 
for  such  predictions.  To  achieve  this  it  is  necessary  to 
identify  the  nonscaling  effects.  One  possibility  is  that 
the  main  source  for  nonscaling  is  the  attenuation  of  the 
blast  wave  due  to  small-scale  roughness.  An  attempt  to  take 
this  factor  into  account  was  presented  in  chapter  6.  However, 
the  time  frame  for  the  project  did  not  allow  an  exhaustive 
analysis,  so  no  firm  conclusion  has  been  reached. 

It  is  believed,  however,  that  the  experimental  data  contained 
in  this  report,  together  with  other  available  information, 
will  make  further  progress  possible.  It  is  therefore  recom- 
mended that  sufficient  resources  should  be  allocated  for  a 
more  detailed  analysis.  The  purpose  would  then  be  prove  if 
the  hypothesis  and  correction  method  suggested  in  chapter  6 
is  correct,  and  to  examine  if  simple  considerations  about 
energy  and  impulse  of  the  blast  wave  can  lead  to  useful  con- 
clusions. 
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Figure  1 The  exper imentc.1  upset,  scale  1:100 
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igure  2 . 1 Sketch  of  tube  with  charge 


51 
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Figure  3.1  CFG  1 in  principle 


Figure  3.2  CFG  2 in  principle 
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Figure  3.3  CFG 
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Definition  of  wall  roughness:  rL 
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Figure  3.5  Definition  of  wall  roughness 
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Figure  3.7  Mean  values  of  front  pressure  versus 
distance  in  three  scales 
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Figure  3.8a  Deviation  from  scaling  with  respect  to  dimension 
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Figure  3.8b  Deviation  from  scaling  with  respect  to  dimension 
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Figure  3 . 8c  Deviation  from  scaling  v/ith  respect  to 


Figure  3.8d  Deviation  from  s caling  with  respect  to 
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igure  5.1a  P0  from  curve  fitting  as  function  of  charge  gruup.  Logarithmic  scale 
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Figure  5.2a  C2  from  curve  fitting  as  function  of  charge  grou;  . Logarithmic  scale 
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igure  5,2b  C2  from  curve  fitting  as  function  of  charge  group.  Logarithmic  scale 
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F igur e 5.3  Difference  in  pressure-time  profile  between 
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Figure  7 . 1 Attempt  on  extrapolation  from  small  scales  to  large  scales 
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q (kgra  3) 

80 

160 

24  0 

320  ! 

r 

j 480 

l 

640 

D=0 . 05m 

10 

20 

30 

40 

60 

80 

Q 

D=0 ,10m 

80 

160 

240 

320 

480 

640 

(g  TNT) 

D— 0 . 15m 

270 

540 

810 

1080 

1620 

2160 

q7v“ 

- *3 

all  D 

12.5 

25.0 

37 . 5 

50.0 

75.0 

100.0 

(kgn 

Table  3.1  Charges  In  the  main  series  of  experiments 


D = 0.15m 

q (kgm  3) 

17.8 

23.7 

29.6 

47.4 

71.1 

Q 

CFG  1 

60 

80 

100 

160 

240 

(g  TNT) 

CFG  2 

80 

100 

160 

240 

Q/V 
kgm  3 

2.78 

3.70 

4.63 

7.41 

11.11 

Table  3.2  Shots  with  small  charges 


e 

50% 

20% 

10% 

5% 

1 

1 2% 

! 

r ' — ] 
1% 

0.2% 

n 

0.674 

1.282 

1.645 

1.960 

2,326 

2,576 

3.090 

Table  4.1  Test  level  and  percentage  points  of 


Student ' s T-distribution 


CFG! 

CFG  2 

rL 

0 

0.03 

0.06 

0 

0.0  3 

0.06 

a:  f 

16  4.4 

160.2 

218.0 

210.0 

97.0 

165.4  | 

1 

T 



5.106 

3.0  46 

3.417 

6 o 687 

2*526 

3.725  1 

Table  4.2a  Values  of  a i and  T in  build-up  region 


CFG1 

CFG  2 

rL 

0 

0.03 

0.06 

0 

0.03 

0.06 

ttj'xlO^ 

24 . 850 

0.770 

-4.405 

-47.500 

-69 . 500 

-71.000 

T 

3.358 

0.091 

0.565 

3.084 

3.971 

4.46  5 

Table  4.2b  Values  of ai1  and  T in  build-up  region 


f'FGl 

1 

CFG  2 

rL 

0 

0.03 

0.06 

0 

0.03 

0.06 

ai  • 

32.8 

18.1 

71.4 

64.2 

32.6 

50.0 

T 

4.530 

3.591 

15.000 

9.093 

6.342 

12.755 

Table  4.3a  Values  of  ax  and  T in  attenuation  region 


CFG1 

CFG  2 

rL 

0 

0.03 

0.06 

0 

— 

0.03 

0.06 

Cii  'xlO3 

0.655 

-2.650 

-7.400 

-9 .600 

-15.650 

-15.650 

T 

0.712 

1.4  56 

2.966 

5*647 

3.525 

4.369 

Table  4.3b  Values  of  gf  and  T in  attenuation  region 
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CFG  1 

CFG  2 

rL 

0 

0.03 

0.06 

i 

0 

0.03 

0.06 

a 2 x 1 0 2 

-2.688 

-2.400 

-]  .338 

-2.788 

-5.588 

-4.900 

T 

3.258 

3.975 

2.352 

2.411 

6.507 

7.582 

a 3 x 1 0 1 

-1.080 

-4.530 

-7.850 

-2.600 

-1.670 

-4.420 

T 

1.306 

7.342 

13.534 

2.031 

1.740 

6.17  3 

Table  4.4a  Values  of  a2 , a 3 and  corresponding  T-values 
in  attenuation  region 


CFG  1 

| CFG  2 

1 

rL  . 

0 

0.03 

0.06 

0 

0.03 

0.06 

0(2  X 10 6 

-2.063 

-23.094 

-49.375 

-10.594 

-29.750 

-50.313 

T 

1.654 

9.146 

14.234 

3.498 

3.437 

7.630 

a 3 x 104 

1.243 

5.525 

10.325 

2.725 

7.050 

11.075 

T 

7.565 

16.617 

22.692 

6.566 

6.104 

12.621 

Table  4.4b  Values  of  a;1  , g?'  and  corresponding  T-values 
in  att enuation  region 
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q = 80 

kgm  3 

CFG 

1 

CFG 

2 

D(m) 

rL 

Cl 

C2 

PD(bar) 

R 

Cl 

C2 

°o  (bar ) 

R 

0 

21.22 

0,61 

56.4 

0 . 91 

15.13 

0.77 

52.4 

0.94 

0.05 

0.03 

6.67 

1.26 

51.1 

0.96 

6.29 

1.05 

33,9 

0.96 

0.06 

2.69 

1.65 

38.4 

0.96 

3.64 

1.23 

26.5 

0.97 

0 

18.91 

0.75 

63,2 

0.96 

17.07 

0.68 

50.8 

0.95 

0.10 

0.03 

6.06 

1.18 

40.5 

0.96 

5.13 

1.09 

29.7 

0.83 

0.06 

4.71 

1.36 

42.3 

0.96 

3.59 

1.33 

30.3 

0.90 

0 : 

21.07 

0.75 

70.1 

0.97 

17.91 

0.81 

66.3 

0.95 

0.15 

0.03 

10.86 

0.86 

43.5 

0.96 

7.25 

1,08 

41.1 

0.95 

0.06 

4.28 

1.51 

48.6 

0.98 

4 .94 

1.15 

31.6 

0.97 

Table  5,1a  Results  from  curve  fitting,  q - 80  kg/m ~ 


q = 160kgm“3 

CFG 

1 

CFG 

2 

D(m) 

rL 

Cl 

C?  : 

?o(bai- ) 

R 

C1 

C2 

:30  (bar) 

R 

0 

18.34 

0 .63 

78.0 

0.77 

16.24 

0.60 

65.0 

0.91 

0.05 

0.03 

5.24 

1.07 

60.9 

0.96 

3.51 

1.13 

47.4 

0.97 

0.06 

1.77 

1.51 

56.9 

0 .96 

2.23 

1.23 

38.2 

0 .98 

0 

19.03 

0.65 

84.7 

0.96 

16.99 

0.58 

64.6 

0 . 89 

0.10 

0.03 

8.97 

0.78 

53.5 

0.94 

4.73 

1.00 

46.7 

0.94 

0.06 

2.25 

1.44 

62.4 

0.97 

3.28 

1.09 

40.6 

0.90 

0 

24.60 

0.57 

90.7 

0.94 

11.41 

0.89 

88.0 

0.89 

0.15 

0.03 

6.26 

0.97 

58.4 

0.97 

6.00 

0.96 

55,3 

0.97 

0.06 

2.39 

1.42 

62.2 

0.97 

2.64 

1.29 

51.5 

0,98 

Table  5.1b  Results  from  curve  fitting,  ct  = 160  kg/m3 
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q = 240kgm  ^ 

CFG 

1 

CFG 

2 

D (m) 

rL 

Cl 

C2 

?o«bar> 

R 

Cl 

C-2 

?o (bar) 

R 

0 

26.93 

0.44 

88.1 

0.87 

20 . 38 

0.49 

77.7 

0.93 

0.05 

0.03 

5.54 

0.93 

69.3 

0.93 

3.82 

1.00 

57.1 

0.96 

0.06 

0.81 

1.56  ' 

55.0 

0.89 

1.05 

1.38 

44.6 

0.93 

0 

29.72 

0 . 36 

7 8.5 

0 . 90 

19.16 

0 . 50 

74.7 

0.91 

0.10 

0.03 

5.12 

0.93 

63.1 

0.94 

3.47 

1.01 

53.2 

0.89 

0.06 

0.93 

1.59 

69.6 

0.87 

1.94 

1.21 

51.4 

0.95 

0 

38.98 

0.27 

82 . 0 

0.77 

24.88 

0.43 

79.7 

0.95 

0.15 

0.03 

6.25 

0 . 88 

68.3 

0.95 

5.02 

0.93 

62.0 

0.95 

0.06 

2.42 

1.23 

> ■ ■ 

67.5 

0.93 

1.83 

1.26 

55.8 

0.96 

Table  5.1c  Results  from  curve  fitting,  q = 240  kg/m^ 


q = 320kgm  ^ 

CFG 

1 

CFG 

2 

D(m) 

rL 

Cl 

c2  : 

P0(bar } 

R 

Cl 

c2  ' 

?Q (bar) 

R 

0 

26.53 

0.42 

92.5 

0.87 

22.69 

0.44 

84.8 

0.94 

0.05 

0.03 

7.05 

0.75 

67.4 

0.92 

4 . 54 

0.85 

58.2 

0.94 

0.06 

0.97 

1.38 

60.6 

0.91 

0.94 

1.36 

55.0 

0.97 

0 

28.51 

0.35 

81.0 

0.65 

23.22 

0.41 

78.3 

0.88 

0.10 

0.03 

6.14 

0.75 

5 7.7 

0 . 86 

3.67 

0.90 

54.5 

0.90 

0.06 

2.05 

1.17 

67.1 

0.90 

1.97 

1.11 

55.9 

. 0.97  . 

0 

43.45 

0.21 

82.1 

0.75 

23.29 

0.46 

91.8 

0.95 

0.15 

0.03 

8.32 

0.69 

65.3 

0.96 

4.40 

0 . 89 

63.7 

0.94 

0.06 

2.59 

1.10 

70.4 

0.91 

1.84 

57.1 

= 320  kg/m^ 


Table  5. Id  Results  from  curve  fitting,  q 
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q = 4 30kgrn  ^ 

CFG 

i 

CFG 

2 

D (m) 

rL 

Cl 

C2 

P0(ba.r) 

R 

Cl 

C2 

?0 (bar) 

R 

0 

18.79 

0.51 

. 105.2 

0.88 

23.38 

0.39 

88.2 

0.78 

0.05 

0.03 

5.58 

0.73 

67.7 

0.93 

3.21 

0.90 

69.6 

0.95 

0.06 

1.57 

1.0  8 ' 

61.6 

0.90 

] .12 

1.16 

57.4 

0.97 

0 

26.22 

0.38 

94.1 

0.91 

21.65 

0.39 

81.9 

0.93 

0.10 

0.03 

4.68 

0.78 

67.5 

0.85 

3.67 

0.90 

73.6 

0.96 

0.06 

2.40 

0.98 

67.0 

0.95 

1.70 

1.03 

56.5 

0.96 

0 

37.83 

0.27 

93.9 

0.87 

24.01 

0.42 

99.2 

0.90 

0.15 

0.03 

7.63 

0.67 

75.1 

0.96 

4.28 

0.80 

65.2 

0.38 

0.06 

3.80 

0.85 

74.3 

0.90 

1.53 

1.10 

64.7 

0.94 

Table  5.1e 


Results  _frorn  curve  fitting 


g = 4 80 


kg/m  ^ 


q = 640kgm  ^ 

CFG 

1 

CFG 

2 

D (m) 

rL 

Cl 

c2  ; 

Po(bar) 

R 

Cl 

c2 

?o (bar) 

R 

0 

15.54 

0.54 

115.7 

0.94 

22.61 

0.37 

89.3 

0.70 

j 

0.05 

0.03 

4.73 

0.75 

75.0 

0.96 

3.41 

0.80 

65.6 

0.79 

:0 . 06 

1.17 

1.12 

71.5 

0.96 

1.13 

1.06 

56.5 

0.91 

0 

25.44 

0.35 

92.5 

0.90 

0.10 

0.03 

5.23 

0.72 

74.9 

0.97 

0.06 

1.32 

0.98 

67.6 

0.94 

0 

27.28 

0.36 

104.7 

0.87 

0.15 

0.03 

8.84 

0.57 

71.1 

0.91 

0.06 

4.12 

0.76 

69.1 

0.90 

Table  5. If 


Results  from  curve  fitting , 


g = 640 


kg/in^ 
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CFG  1 Po  “ ac2 


D (m) 

rl'i 

a 

b 

R 

p0 (q~8o ) 

0 

13.8 

0.33 

r 

0.99 

59.0 

0 . 05 

0.03 

25.9 

0.17 

0.93 

53.3 

0.00 

13.5 

0.26 

0.93 

41.7 

0 

31.8 

0.17 

0.89 

66.9 

0.10 

0.03 

12.9 

0.27 

0.96 

42.5 

0.06 

20.0 

0.20 

0.81 

48.7 

0 

36.8 

0.15 

0.84 

72.3 

0.15 

0.03 

16.4 

0.24 

0.92 

47.2 

0.06 

23.8 

0.18 

0.89 

52.4 

Table  5.2a  Dependence  of  Pp  on  q 


CFG  2 

Po  = 

= aqb 

D (m) 

rL 

a 

b 

R 

Po (q=80) 

0 

16.8 

0.27 

0.97 

54.6 

0.05 

0.03 

8.5 

0.33 

0.96 

36.4 

0.06 

5.3 

0.39 

0.96 

28.5 

0 

13.8 

0.30 

0.98 

51.6 

0.10 

0.03 

3.4 

0.50 

0.98 

30.5 

0.06 

6.1 

0.37 

0.97 

31.5 

■ 

0 

28.2 

0.20 

0.90 

68.8 

0.15 

0.03 

13.8 

0,26 

0.95 

43.6 

0.06 

6.6 

0.38 

0.94 

1 34.8 

Table  5.2b  Dependence  of  Pp  on  q 
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C2  = a + bD 


CFG  1 

' ; 

f 

j 

L - 0 

rL  = 0 

. 03 

*1 

, = 0.06 

q (kgm~3) 

a 

b 

R 

a 

b 

R 

a 

b 

R 

80 

0.56 

1.4 

0.87 

1.50 

-4.0 

0.94 

1.65  1 

-1.4 

0.48 

160 

0.68 

-0.6 

0.72 

1.04 

-1.0 

0.35 

1.55 

-0.9 

0.95 

240 

0.53 

-1.7 

1.00 

0.96 

-0.5 

0.87 

1.79 

-3.3 

0.82 

320 

0.54 

-2.1 

0.98 

0.79 

-0.6 

0.87 

1.50 

-2.8 

0.96 

480 

0.63 

-2.4 

1.00 

0.79 

-0.6 

0.55 

1.20 

-2.3 

0 . 9 9 

640 

0.60 

1.,.-  . ■ „ 

-1.8 

0.84 

0.86 

-1.8 

0.93 

1.31 

-3.6 

0.99 

Table  5.3a  Dependence  of  C2  on  D 


C 2 = a f bD 


CFG  2 

1 

i 

t"1 

!l 

0 

rL 

= 0.03 

rL  = 0.06 

q (kgm-3) 

a 

b 

R 

a 

b 

R 

a 

b 

R 

80 

0.71 

0.4 

0.30 

1.04 

0.3 

0.72 

1.32 

-0.8 

0.45 

160 

0.40 

2.9 

0.84 

1.20 

-1.7 

0.95 

1.14 

0.6 

0.30 

240 

0.53 

-0.6 

0.79 

1.05 

-0.7 

0.80 

1.40 

-1.2 

0.69 

320 

0.42 

0.2 

0.40 

0,84 

0.4 

0.75 

1.42 

i-2,1 

0.78 

480 

0.37 

0.3 

0,87 

0.97 

-1.0 

0.87 

1.16 

-0.6 

0.46 



Table  5 » 3b  Dependence  of  C2  on  D 
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b 

= b],q  + 

b2 

CFG 

1 

CFG 

2 

rL 

b q • 1 0 3 

b2 

R 

bl* 

103 

b2 

R 

3 

-4.9 

0. 

36 

0 

.72 

-2 

.7 

i. 

34 

0 

.32 

0.03 

2.2 

-2 . 

13 

0 

.35 

-0 

.8 

-o. 

33 

0 

.14 

0.06 

-3.6 

-1. 

24 

0 

.70 

< 

"I 

I 

* 8 

-o. 

37 

0 

.28 

Table  5.4  a and  b Interaction  between  charge  group  and  scale 


„ b 

C2  = aq 


CFG  1 

CFG  2 

D (m) 

rL 

a 

b 

R 

a 

b 

R 

0 

0.91 

-0.10 

0.46 

" 3.76 

-0.37 

0.99 

0.05 

0.03 

4.35 

-0.28 

0.95 

L'  2.13 

-0.15 

0.84 

0.06 

4.40 

-0.21 

0.89  1 

1.67 

-0.06 

0.42 

0 

4.15 

-0.40 

0.87  1 

1.38 

-0.18 

0.62 

0.10 

0.03 

2.60 

-0.21 

0.82 

j 1.78 

-0.11 

0.93 

0.06 

3.73 

-0.20 

0.74 

2.19 

-0.12 

0.82 

0 

5.05 

-0.47 

0.73 

; 6.39 

-0.4  5 

0.84 

0.15 

0.03 

2.60 

-0.22 

0.82 

' 2.16 

-0.16 

0.9  8 

0.06 

7.70 

-0,35 

0.95 

1.41 

-0 . 03  j 

0.33 

Table  5.5  a and  b Dependence  of  C2  on  a 
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= a + brL 


CFG  1 

D = 

0.05 

m 

D = 

o 

i — i 

O 

ii 

m 

D 

= 0.15 

m 

q (kgm-3) 

a 

b 

R 

a 

b 

R 

a 

b 

R 

80 

0.65 

17.3 

0.99 

0.79 

10.2 

0.97 

0.66 

12.7 

0.93 

160 

0 . 63 

14.7 

1.00 

0.56 

13  ..2 

i 

0.93 

1 

0.56 

14.2 

1.00  | 

240 

0.42 

18.7 

0.99 

0 . 35 

20.5 

! l.oo 

0.31 

16.0 

0.99  1 

320 

0.37 

16.0 

0.98 

0 . 35 

13.7 

1.00 

0 .22 

14.8 

1.00  | 

480 

0.49 

9.5 

0.99 

0 .41 

10.0 

0,90 

0 . 31 

9.7 

0.97  | 

640 

0.51  ] 

9.7 

0.99 

0.37 

10.5 

0.99 

0.36 

6.7 

1.00  j 

Table  5.6a  Dependence  of  C2  on  rT, 


C2  = a + brL 


CFG  2 

D = 

0.05 

D 

•=  0.10 

D = 

= 0.15 

q (kgm-3) 

a 

b 

R 

a 

b 

R 

a 

b 

R 

80 

0.79 

7.7 

0.99 

0,71 

10 .8 

0 . S 9 

0.84 

5.7 

0.95 

16  0 

0.67 

10.5 

0.93 

0.64 

8.5 

0.94 

0.84 

6.7 

0.9  4 

240 

0.51 

14.8 

0.99 

0.55 

11.8 

0.97 

0.46 

13.8 

0.99 

320 

0.42 

15.3 

1.00 

0.46 

11.7 

0.97 

0.49 

11.5 

0.99 

480 

0.43 

12.8 

0.98 

0.45 

10 .7  1 

0.94 

0.43 

11.3 

1.00 

640 

0.40 

11.5 

0.99 

1 

J 

Table  5.6b  Dependence  of  C 2 on  x~l 
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CFG  1 


i = 

1 

2 

3 

4 

5 

6 

7 

1.09 

2.0 

27.6 

-0.094 

-20 

-0.55 

-2.2 

F-i 

0.32 

0.90 

0.21 

0.26 

0.15 

0.34 

Tahle  5.7a  Results  from  linear  regression  on  C2 


CFG  2 


i = 

1 

2 

3 

4 

5 

6 

n 

/ 

»± 

1.47 

2.75 

-6.16 

-0.175 

-24 

-0.42 

2.7 

Ri 

1 

0.23 

0.23 

0.18 

0.16 

0.05 

0.19 

Table  5 . 7b  Results  from  linear  regression  on  C? 
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CFG  1 

rL  = 0 

D (m) 

curve 

A 

;i 

CO 

o 

q-160 

q-240 

q=320 

q=480 

q=64  0 

i 

2.17 

1.44 

1.25 

1,53 

1.20 

1.02 

0.05 

ii 

3.04 

1.09 

1.02 

2.62 

2.41 

2.01 

i 

1.61 

2.81 

1.46 

2.41 

7.14 

1.50 

0.10 

ii 

2.83 

78.56 

2.80 

2.23 

4.50 

1.73 

T 

1.25 

1.12 

1.08 

2.09 

1.99 

3.82 

0.15 

II 

1.10 

2.02 

1.06 

1.95 

3.74 

2.19 

i 

Table  5.8  Comparison  between  curveforms 


D = 

0.15 

° 

11 

u 

CFG 

Q(g) 

q (kgm 

Cl 

C2 

Po (bar) 

R 

60 

17.8 

22.54 

0.96 

24.9 

0.99 

80 

23.7 

21.29 

0.95 

31.0 

0.99 

1 

100 

29.6 

20.68 

0.92 

| 

1 36.4 
| 

0.98 

160 

47.4 

20.87 

0.77 

4 8.1 

0.98 

240 

71.1 

21.00 

0.72 

61.4 

0.98 

80 

23.7 

19.59 

0.82 

27.0 

0.75 

2 

100 

29.6 

21.90 

0.89 

37.9 

0.9  4 

160 

47.4 

18.91 

0 . 84 

47.2 

0.98 

240 

71.1 

14.70 

0.83 

51,0 

0.95 

Table  5.9  Results  from  curve  fitting 
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Table  6.1  Shock  velocities 
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Table  6.2  Values  of  Zliliim  and  LOI-A°0> 

P ' (1:100)  P' (2:100) 
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SHOT  2 

Q = 100  kg  g = 

__  

= 5.4  kgm 

n 

7.1 

10.9 

14.7 

18.4 

P ''bar) 

1 

12.6 

I 

8.4 

! 

5.4 

5.2 

< 

SHOT  6 

Q = 

300  kg 

q = 16.1  kgm 

-3 

n 

7.1 

10.9  j 

: 14.7 

18.4 

26.1 

29.9 

P (bar 

26.4 

15.9  j 

1 

12.7 

7.1 

6.7 

2.7 

SHOT  8 

Q = 

1000  kg 

q = 53.7  kgm 

-3 

n 

8.3 

12.0 

15.8 

19.6 

27.3 

31.0 

P (bar) 

59.5 

40.7 

19.0 

15.7 

13.7 

6.0 

SHOT  9 

Q » 1000  kg 

q = 53.7 

kgm  ^ 

n 

6.7 

14.2 

H 

CO 

9 

O 

25.7 

29.5 

P (bar) 

85.9 

20.7  j 

! 18.5 

12.5 

4.4 

Table  7.1a  Results from  Raufoss  exp e rime n t s 


CFG  1 

D = 0*15  r l = 0 

17.8 

23.7 

29.6 

47.4 

71.1 

4 

78.8 

86.0 

93.7 

108.8 

113.1 

8 

47.6 

58.4 

63.3 

79.1 

82.3 

12 

29.1 

37.3 

42.0 

53.2 

64.6 

16 

25.6 

30.1 

35.6 

46.5 

59.7 

20 

19.7 

25.3 

29.3 

39.8 

51.2 

24 

16.0 

20.9 

25.4 

34.6 

45.5 

28 

15.2 

18.8 

22.5 

33.3 

43.6 

32 

12.9 

16.5 

20.2 

29.8 

39.8 

Table  7.1b  Some  results  from  model  tests,  small  charges 
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CFG  1 

q : 

= 80  kgm 

-3 

q = 

160  kgm 

-3 

r. 

\D  (m) 

rL  " \ 

0.05 

0.10 

0.15 

0.05 

0.10 

r~ 

0.15 

0 

87.4 

95.5 

83.4 

98.6 

67.2 

88.2 

4 

0.03 

S0.1 

62.3 

82.1 

83.7 

37.3 

77.7 

0.06 

42.5 

53 . 7 

80.7 

75.6 

64 . 4 

97.8 

0 

74.0 

71.4 

77.3 

87.9 

} 

82.3 

99.3 

8 

0.03 

68.0 

47.6 

72.2 

i 81.7 

66.1 

79.9 

0.06 

45.3 

46.9 

60.9 

! 67.6 

58.1 

72.6 

r 

0 

6 6.6 

64.3 

81.3 

83.1 

79.6 

98.7 

12 

0.03 

50.9 

40.3 

56.5 

61.1 

51.3 

6 3.9 

0.06 

46.6 

35.9 

49.2 

60.8 

44.5 

53.7 

0 

53.7 

61.4 

63.0 

72.4 

86.3 

89.2 

16 

0.03 

48.7 

34.6 

39.9 

55.6 

50.1 

57.3 

0.06 

32.7 

38.0 

45.3 

45.5 

55.9 

53.8 

0 

49.9 

55.3 

64.6 

74.0 

76.8 

89.5 

20 

0.03 

3 8.1 

31.9 

34.8 

48.3 

44.5 

45.0 

0.06 

30.9 

33.6 

32.1 

43.3 

42.7 

4 2.7 

0 

46.5 

48.4 

53.0 

54.0 

66.4 

67.7 

24 

0.03 

33.2 

28.0 

34.4 

39.7 

38.8 

38.0 

0.06 

18.1 

27.6 

29.7 

36.0 

39.8 

40 . 9 

0 

39.9 

44.2 

50.2 

61.5 

56.9 

64.1 

28  ; 

0.03 

38.2 

23.0 

28.2 

37.8 

34.3 

36.5 

0.06 

17.3 

20 ,1 

23.2 

29.4 

30.7 

31.0 

0 

36.3 

37.3 

39.8 

53.1 

50.2 

58.0 

32 

0.03 

19.8 

20 . 2 

24.4 

30.3 

36.9 

32.5 

0.06 

12.5 

16.6 

18.4 

20.7 

24.5 

26.3 

Table  7.1c  Some  results  from  model  tests,  main  series 
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(m) 

q ( kgm” 3 ) ' ■ . 

0.05 

0.10 

0.15 

2.65 (Raufoss) 

5.4 

2.19 

2.17 

2.14 

1.00 

16.1 

1.90 

1.85 

1.79 

1.40 

53.7 

1.60 

1.51 

1.42 

1.75  (1.57) 

Table  7.2  Values  of  C?  obtained  from  the  Raufoss 
experiments  and  from  the  model  tests 


Q = 10  kg  q = 0.16  kgm  ^ 

n 

2.25 

5.00 

8.75 

14.50 

P (bar) 

1.9 

1.4 

1.0 

0.4 

Q = 57  kg  q = 0.89  kgm  ^ 

h 

2.0 

5.0 

10.0 

14.5 

20.0 

P (bar) 

6.3 

3.7 

2.0 

1.8 

1.3 

Q=166kg  q = 2.6  kgm 


n 

5.0 

7.5 

12.5 

22.5 

P (bar) 

9.3 

7.5 

4.0 

2.3 

~ 3 " 

Q =500  kg  q = 7.8  kgm 

n 

3.5 

5.0 

6.5 

10.0 

20.0 

27.5 

P (bar) 

34.0 

21.0 

13.6 

6.9 

3.8 

2.6 

Table  7.3  Results  from  Lista  experiments 
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(m) 

q(kgm“3p\v^ 

0.05 

0.10 

0.15 

4.0  (Lista) 

0.16 

3.26 

3.35 

3.44 

0.78 

0.89 

2.81 

2.86 

2.90 

0.68 

2.6 

2.54 

2.56 

2.57 

0.97 

7.8 

2.26 

2.25 

2.23 

1.24 

Table  7.4  Values  of  C2  obtained  from  the  Jjista 


experiments  and  the  model  tests 


